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Abstract 


Existing interior noise reduction techniques for aircraft 
fuselages perform reasonably well at higher frequencies, but are 
inadequate at lower frequencies, particularly with respect to the 
low blade passage harmonics with high forcing levels found in 
propeller aircraft. This research focuses on a noise control method 
which considers aircraft fuselages lined with panels alternately 
tuned to frequencies above and below the frequency that must be 
attenuated. Adjacent panels would oscillate at equal amplitude, to 
give equal source strength, but with opposite phase. Provided 
these adjacent panels are acoustically compact, the resulting 
cancellation causes the interior acoustic modes to become cutoff, 
and therefore be non -propagating and evanescent. This interior 
noise reduction method, called Alternate Resonance Tuning (ART), 
is described in this thesis both theoretically and experimentally. 

Problems presented herein deal with tuning single paneled 
wall structures for optimum noise reduction using the ART tuning 
methodology, and three theoretical problems are analyzed. The 
first analysis is a three dimensional, full acoustic solution for tuning 
a panel wall composed of repeating sections with four different 
panel tunings within that section, where the panels are modeled as 
idealized spring-mass-damper systems. The second analysis is a 
two dimensional, full acoustic solution for a panel geometry 
influenced by the effect of a propagating external pressure field 
such as that which might be associated with propeller passage by a 
fuselage. To reduce the analysis complexity, idealized spring-mass- 
damper panels are again employed. The final theoretical analysis 
presents the general four panel problem with real panel sections, 
where the effect of higher structural modes is discussed. Results 
from an experimental program highlight real applications of the 
ART concept and show the effectiveness of the tuning on real 
structures. 
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Chapter 1 


The General Noise Problem 
and 

Alternate Resonance Tuning 


Section l.i: Introdttcti.an 

Existing noise reduction techniques for paneled structures 
perform reasonably well at higher frequencies, but are less 
effective at lower frequencies, particularly if the low frequency 
noise problem has high forcing levels such as those found in 
helicopter or propeller aircraft. This research effort focuses on a 
method which considers enclosures lined with panels alternately 
tuned to frequencies above and below the frequency that must be 
attenuated. Adjacent panels would oscillate at equal amplitude, to 
give equal source strength, but with opposite phase. If these 
adjacent panels are acoustically compact, the resulting 
cancellation causes the interior acoustic modes to become cutoff, 
and therefore be non -propagating and evanescent. This noise 
reduction method, called Alternate Resonance Tuning (ART), has 
potential application to reducing interior noise in a variety of 
structures. As mentioned previously, one particularly appropriate 
application is in aircraft fuselages. With some minor 
modifications, the method might also be applied to a variety of 
noise problems in other areas, such as road and wind noise in 
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automobiles, flow and fan noise in heating and air conditioning 
ducts, and noisy discrete frequency components in household 
appliances. 

The ART technique is a procedure intended to reduce low 
frequency noise within an enclosure. A paneled wall could be 
constructed of, or lined with, a series of special panels which 
would allow the designer to control the wavenumber spectrum of 
the wall motion, thus controlling the interior sound field. By 
judicious tuning of the structural response of individual panels, 
wavelengths inside the enclosure can be reduced to the order of 
the panel size, thus causing low frequency interior acoustic modes 
to be cutoff provided these panels are sufficiently small. By 
cutting off the acoustic modes in this manner, a significant 
reduction of interior noise at lower frequencies should be 
achieved. 

Current aircraft fuselage noise control treatments have 
already demonstrated that the mass and stiffness of individual 
panels can be altered. This research effort demonstrates that 
panel resonance frequencies can be manipulated to achieve the 
ART effect. Application of this concept might involve the 
modification of existing structural panels or development of a new 
design for enclosure interior trim panels. Although complete 
acoustic cutoff will not be achievable in practice, an approximate 
cancellation should still substantially reduce the interior noise 
levels at the particular frequency of interest. It is important to 
note that the ART method utilizes the flexibility and dynamic 
behavior often found in these structures to good advantage, 
although these properties are not normally beneficial in noise 
control. 

Section 1.2; Review of the Literature 

A survey of the current state-of-the-art in noise 
measurements and noise reduction techniques associated with 
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paneled structures provides insight for the problems found in 
these structures. The following is not intended as a complete 
literature review but rather a survey of work related to the 
definition of the noise problem in paneled structures (especially 
aircraft applications) and existing techniques for noise reduction in 
those structures. Low frequency sound transmission of propeller 
near-field pressure, radiated noise, and structural vibration • 
through the aircraft fuselage are important sources of interior 
noise for propeller aircraft. Test results of interior noise 
measurements on actual aircraft, such as those conducted by Kuntz 
and Prydz 1 and Wilby and Wilby 2 , confirm the presence of high 
noise levels outside and inside the fuselage. Figure 1.2.1 
(reproduced from Kuntz and Prydz) shows sound pressure level 
spectra measured outside the fuselage behind the propeller plane 
with and without the propfan in place. Figure 1.2.2 (also 
reproduced from Kuntz and Prydz) shows the resulting interior 
cabin spectra measured with and without the propfan. There 
remains a definite need to reduce the interior noise due to this 
discrete, low frequency noise source. The ongoing development of 
advanced turboprop aircraft technology for commercial use 
(spurred forward by fuel savings claims of 35% for a 
contrarotating propfan over a conventional turbofan 3 ) underscores 
the need for a solution to this noise problem. The first few 
propeller blade passage harmonics are difficult to reduce because 
of their low frequency and high forcing levels. Low frequency 
acoustic transmission is particularly difficult to block or absorb 
given the weight and wall thickness constraints imposed by 
fuselage construction techniques. Fuselage panels are necessarily 
thin and flexible and are often lightly damped, with resonance 
frequencies in the same range as the blade passage harmonics. 
Helicopters suffer from similar discrete low frequency noise 
problems, where the blade passage frequencies are even lower 
(and hence, harder to control) than that found in propeller driven 
aircraft, as shown in Figure 1.2.3 (reproduced from Brentner). 4 
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Frequency (Hz) 


Figure 1.2.1: Exterior fuselage sound pressure spectra measured 63.5 cm 

behind the propeller plane with and without the propfan on a Gulfstream II 
aircraft (reproduced from Kuntz and Prydz). 



Frequency (Hz) 


Figure 1.2.2: Cabin sound pressure level spectra measured in the prop plane 

with and without the propfan (reproduced from Kuntz and Prydz). 
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Figure 1.2.3: Experimental data for the first 40 harmonics of the blade 

passage frequency taken from 1/4 scale model tests, including a comparison 
to predicted levels by the noise prediction program WOPWOP (reproduced 
from Brentner). 

Although not mentioned as often as discrete noise sources, 
broadband noise can also play a role in fuselage interior sound 
levels at low speed, low power operating conditions, as reported by 
Knowles . 5 

An additional and extensive review of the subject of 
propeller-driven aircraft interior noise is given by Mixson and 
Powell . 6 The more specific topic of noise transmission through 
aircraft panels is reviewed by Vaicaitis, Grosveld, and Mixson . 7 A 
current methodology to predict the performance of acoustic 
treatments installed in an actual aircraft fuselage is given by 
Heitman and Mixson . 8 These papers show the complexity of the 
interior noise problem and demonstrate the degree of success of 
current acoustic treatment methods. In the detailed design of an 
interior noise treatment, as described by Vaicaitis and Mixson^, 
panels have had alterations in mass, stiffness, and damping, and 
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fuselage walls have been filled with multi-layered absorptive 
materials. Figure 1.2.4 (reproduced from Vaicaitis and Mixson) 
displays the breadth and complexity of such acoustic treatments. 

It has been found most effective to utilize different treatment 
combinations depending on the specific panel or panel region 
involved. Considerable noise reduction has been achieved over 
most of the frequency spectrum by careful application of these 
standard methods. However, additional improvement is still 
required at the lowest harmonics of blade passage frequency, 
where even sophisticated multi-layer wall treatments are not 
adequate. Figure 1.2.5 (also reproduced from Vaicaitis and 
Mixson) shows a comparison of measured noise levels in an 
untreated aircraft and a prediction of noise levels after acoustic 
treatment has been applied. Note that little improvement is 
realized below about 200 Hertz. 

Beyond the implementation of these basic noise reduction 
techniques (adding mass, stiffness, structural damping, and 
absorptive material), several other approaches specifically 
applicable to this problem have been previously put forth, or are 
currently being investigated. These include the careful design of 
the periodic structural members of the fuselage and the concept of 
intrinsic structural tuning. 10.1 1.12 Another approach involves the 
use of a large number of sharply tuned Helmholtz resonators 
embedded within the fuselage wall. 13 * 14 Other techniques include 
direct consideration of the propeller as a noise source. One such 
possibility is the reduction of propeller noise by controlling the 
phasing of multiple propellers. 6 Another promising area of current 
research is the use of active noise control. 1 ^* 16 This method seeks 
to achieve cancellation of certain acoustic modes in the cabin 
enclosure by the introduction of additional canceling acoustic 
sources. These sources are intended to actively adjust their 
properties to achieve the desired effect. 

There is a clear need for a new method to reduce discrete, 
low frequency propeller generated interior noise. Such a method 
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Figure 1.2.4: Distribution of acoustic treatments for noise control in an 

aircraft fuselage (reproduced from Vaicaitis and Mixson), 
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Figure 1.2.5: Measurement of noise transmitted into untreated fuselage 

compared with prediction of noise transmitted into treated fuselage 
(reproduced from Vaicaitis and Mixson). 
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must successfully coexist with the noise reduction methods already 
demonstrated to work well at higher frequency, and must also 
function within the constraints imposed by fuselage construction 
techniques. It would be highly desirable to have a method that 
does not entail a significant weight penalty, does not require 
precise tuning, and is relatively robust in terms of variations in 
operating conditions and environment. Alternate Resonance 
Tuning (ART) is a novel approach to low frequency noise reduction 
that appears to satisfy these constraints and still be compatible 
with current noise reduction techniques. This new approach could 
be implemented on existing fuselage structures. Even greater 
effectiveness may be achieved with structures intentionally 
designed to utilize this method, such as future turboprop designs 
or state-of-the-art fuselage trim panels and liners. 

Section 1.3: Alternate Resonance Tuning 

The concept of Alternate Resonance Tuning was developed at 
Duke University by Dr. Donald B. Bliss in 1986.^ Many noise 
control methods depend heavily on the addition of structural 
damping as a means to reduce the motion of surfaces which 
radiate noise. ART relies instead on designing paneled structures 
which use acoustic cutoff as a means to couple poorly to the 
acoustic fields within an enclosure. The method therefore allows 
surfaces to move, but in a motion which is prescribed by the 
designer. Additional damping is not required to improve noise 
reduction performance in a structure provided the structure 
follows the criteria outlined below. 

Consider a low frequency noise source, such as a propeller 
blade passage harmonic, that needs to be suppressed in the 
fuselage interior. The forcing frequency will be fixed by the 
characteristics of the noise source, and the wavelength will 
typically be large compared to most fuselage construction features; 
for example, panel sizes, frame spacings, and so forth. The 
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characteristics of the interior sound field depend on the possible 
modes of motion of the interior wall surfaces. The fuselage walls 
could be constructed using a series of panels. These may be panels 
normally used in fuselage construction, or a special panel structure 
can be introduced on the interior walls. The paneled structure of 
the fuselage walls allows the opportunity to control the 
wavenumber spectrum of the wall motion and thereby control the 
behavior of the interior sound field. Specifically, by judicious 
tuning of the structural response of the individual panels, the 
wavenumber spectrum of the interior wall motion can be reduced 
to the order of the panel size. At low frequency, physically 
reasonable panel sizes correspond to wavelengths for which 
interior acoustic modes are cutoff; i.e., the modes are evanescent 
and nonpropagating. Achieving cutoff will produce a dramatic 
reduction in interior noise levels at the propeller blade passage 
harmonic. 

To achieve cutoff, the panels must be smaller than the 
acoustic wavelength at the frequency in question; namely, the 
adjacent panels must be acoustically compact. Furthermore, the 
motion of adjacent panels must be out of phase and have equal 
and opposite acoustic source strength. For equal size panels, it is 
thus necessary for adjacent panels to execute equal and opposite 
motions. Figure 1.3.1 shows qualitatively the velocity magnitude 
and phase relationships for the lowest modes of two appropriately 
tuned panels. To achieve this condition, the adjacent panels must 
be tuned to have resonance frequencies alternately above and 
below the frequency to be attenuated. For this reason the method 
is called Alternate Resonance Tuning (ART). Although 
complete acoustic cutoff will not be achieved in practice, an 
approximate cancellation will still substantially reduce interior 
noise levels at and around the particular frequency of interest. 

It is important to note that the ART method utilizes the 
flexibility and dynamic behavior of the structure to good 
advantage, although these properties are not normally beneficial in 
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Figure 1.3.1: Amplitude and phase relations for adjacent 

panels for ART tuning. 
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noise control. Whereas traditional noise control methods focus on 
reducing the amplitude of wall motion, the ART method focuses on 
controlling the shape of the wall motion to break the coupling to 
the acoustic interior. Interestingly, the ART method tends to work 
better with less damping, since more perfectly out-of-phase panel 
motions can be achieved. Transmission loss is obtained by 
reflecting the incident sound field, although the individual panels 
are neither massive nor stiff. 


Chapter 2 

The Four Panel Problem 


Section 2.1; Introdiictton 

Many paneled structures are composed of a finite number of 
elements or building blocks which are grouped together into a 
repeating pattern which covers a structure over a desired area. 
Some familiar examples, as noted in Chapter 1, might include an 
aircraft fuselage with panels applied over ribs and stringers, a 
paneled wall on a railroad car or passenger vehicle, a noise 
reducing divider in a building between noisy equipment rooms 
and quiet office areas, and a heating/ventilation duct. The 
theoretical models presented in this chapter consider sound 
transmission through such a wall of panels, as shown in Figure 
2.1.1. A wall of infinite extent composed of idealized spring-mass- 
damper panels is arranged in a periodic pattern. Note that these 
idealized panels are assumed to have a flat, rigid panel element. 
Sound waves are assumed to strike the wall at normal incidence 
and the acoustic transmission is calculated in terms of the dynamic 
properties of the panels. In addition, the analysis considers 
placement of a parallel barrier beyond the panel wall of specified 
acoustic impedance to simulate reflections within an enclosure. 

The panel wall itself is subdivided into identical blocks of panels; 
within each block the analysis allows for four panels with varying 
dynamic properties. By adjusting the panel dimensions 
(wavelengths must be long compared to the panel dimensions) and 
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Figure 2.1.1: Panel wall layout with analysis subsection 

highlighted with thickened line. Note periodicity assumed 
in the wall construction. Panels are labeled in matrix-like 
notation. 

choosing different dynamical properties for the panels, a variety of 
one, two, and three dimensional physical configurations may be 
simulated as shown in Figure 2.1.2. The actual number of degrees 
of freedom depends on the nature of the velocities of the fluid 
close to the surface of the panel array. In Figure 2.1.2, the first 
three dimensional pattern will have a horizontal velocity as fluid 
sloshes back and forth between panels 1 and 2 and panels 3 and 4. 
A vertical velocity is created by similar vibration of panels 1 and 3 
and panels 2 and 4. All panels produce fluid motion in a direction 
perpendicular to the page by virtue of their oscillation. The 
paneled pattern in Figure 2.1.2 repeats indefinitely on the wall, 
and the analysis proceeds with a repeating portion of the four 
panel geometry extracted and placed in a duct, as shown in Figure 
2.1.3. 
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Figure 2.1.2: Arrangement possibilities for 

analysis configuration. 























Section 2.2: Acoustic, Branch Analysis 


The acoustic branch analysis technique is a quasi one- 
dimensional treatment which provides a fairly easy method for 
obtaining approximate solutions for the behavior of a system as 
shown schematically in Figure 2.2.1. Additionally, this method 
provides simple results with which more elaborate analyses may 
be checked. A brief derivation of the branch analysis result will 
introduce terminology, acoustics results, and 
nondimensionalization which will be applied to a variety of 
problems throughout this work. 

With respect to Figures 2.2.1 and 2.2.2, capital subscripts are 
applied to pressures on the exterior side of the paneled structure, 
while small subscripts are applied to pressures on the interior. 

Note that Pi is the uniform normal incident pressure applied to the 
panels, and Pr is the reflected exterior pressure amplitude. 
Similarly, Pt is the pressure transmitted through the structure and 
present on the downstream side, and Pr is the pressure amplitude 
of the reflected wave from the termination impedance Zb- All 
pressures are complex quantities. Other alternative pressure 
variables apply for different analysis purposes; Pa is the net 
effective exterior pressure at point A, while Pb is net effective 
interior pressure at point B on the downstream (right of the panel 
barrier) side. The duct is of length L from panel barrier to 
termination impedance, with a coordinate system located as shown 
in Figure 2.2.2. The analysis will solve for the noise reduction 
which occurs across the panel barrier. This is equivalent to the 
sound pressure level difference in decibels which microphones 
would read at locations Mi and M 2 ; a positive noise reduction is 
indicated by a lower sound pressure level at position M 2 . The 
acoustic loading due to the cutoff modes (inertance or apparent 
mass effects) acting on the panels is neglected in the branch 
analysis, which significantly simplifies the results. Essentially, the 
incident acoustic field is split into separate branches, one for each 
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Figure 2.2.1: Schematic of analysis configuration 

placed in the duct. 

panel as shown in Figure 2.2.2, and these branches are rejoined on 
the other side of the panel barrier. At each branch junction 
pressure and flow continuity conditions are applied. Shape details 
of the individual branches do not affect the results; therefore, a 
coordinate system in the plane of the panel barrier is not required 
for this analysis. 

For the branched panel barrier, note that duct areas on both 
sides on the panel barrier are the same; that is, 

Aa = A® = X Aj (2.2.1) 

i = 1 
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Figure 2.2.2: Equivalent branch analysis schematic. 

where the subscripts denote upstream (A) or downstream (B) 
position with respect to the panel wall barrier. N is the total 
number of branches considered, or 4 for this analysis. A one 
dimensional continuity equation may then be written as 

Vo e iat A a = Uq e i(Bt A B = I Ui e iwt Aj (2.2.2) 

i = 1 

The zero subscript on velocities U represents the one dimensional 
acoustic mode, or the “zero mode”. Let Zmi denote the mechanical 
impedance of the i tlj panel. This mechanical impedance is of the 
form 

Zn.-Ri + ilMiffl + D (2.2.3) 

where Ri is the panel damping constant. Mi is the panel mass, and 
Si is the panel stiffness. Figure 2.2.3 shows a typical 
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Figure 2.2.3: Typical nondimensional impedance 

versus nondimensional frequency with components 
as noted. Low damping yields small but nonzero 
real component of the impedance. 


nondimensional impedance versus frequency plot. Note that the 
magnitude of the nondimensional impedance approaches a 
minimum at the panel nondimensional natural frequency, and that 
the real part of this impedance is frequency independent. 

Invoking the definition of mechanical impedance along with a 
simple force balance on the panel barrier yields 


UiZ m , = Fi = (Pj - Po)Ai (2.2.4) 

Solving Equation (2.2.4) for Ui and inserting into Equation (2.2.2) 
yields 



Taking advantage of the fact that Aa = Ab and 
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(2.2.5) 



,,b_p£ 

u o -~r 
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the following relationship may be obtained. 



( 2 . 2 . 6 ) 


Here Zo B refers to the acoustic impedance looking downstream 
from the panel barrier. A derivation of this acoustic impedance 
Zo B at the point X = -L begins by assuming right and left traveling 
pressure and velocity expressions given by 

p(x,t) = Ae - kx > + Be + kx > (2.2.7a) 

u(x.t) = A* i(ot - kx) _ B e + kx) (2.2.7b) 

v pc pc 

Standard acoustic impedance at x and t is defined as 

Zb = E (M) (2.2.8) 

u(x,t) 

At x = 0, the following classic result is obtained after substitution 
into Equations (2.2.7). 

l = Zb .~ P i (2.2.9) 

A Zb + pc 

Substitution into (2.2.7) for X = -L followed by substitution of 
(2.2.7a) and (2.2.7b) into (2.2.8), coupled with the use of (2.2.9) 
and some manipulation produces 
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( 2 . 2 . 10 ) 



Zb cos kL + ip 0 c sin kL 
p 0 c cos kL + iZb sin kL 


Equation (2.2.6) may be then solved for the ratio of downstream to 
upstream pressures adjacent to the panel barrier in Figure 2.2.2 as 


Pg_ 

Po Aa. + l Al 

2? i= 1 Zn, 


( 2 . 2 . 11 ) 


However, it is more useful to determine the effective pressure at 
the microphone location -Lm’» in this manner, the analysis will yield 
results which are equivalent to a microphone reading in decibels at 
position -Lm- Note that 


M 


p® p® p® 

pA 
M) 


pA pB 
^0 ^0 


( 2 . 2 . 12 ) 


where 


P Zb cos k[L - Lm] + ipoC sin k[L - Lm] ^2 2 13) 

pB Zb cos kL + ip D c sin kL 

Equation (2.2.13) is obtained in a similar manner as was Equation 
(2.2.10). Equation (2.2.12) may then be directly turned into a 
noise reduction in decibels as 


NR = -10 log 


p® 

iJL 

pA 

*0 


-10 log 


P B 

LM 

P B 


(2.2.14) 
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A similar methodology can determine the ratio of transmitted to 
incident pressure as 


(2.2.15) 


(2.2.16) 


Here Zq a is the impedance on the left hand side of the panel 
assembly. Equation (2.2.12) may be applied for any number of 
panels to obtain closed form dimensional results. For example, for 
2 panels and an anechoic termination, 

Ei = Z| + .& (2.2.17) 

Pi Zi + Zi + ? 1 -% L 

PocWH 

For four panels and an anechoic termination. 

Pi _ z» z 2 z 3 -i- z 2 z 3 Z4 + Z 3 Z4 Zi + Z4 Zi z 2 

Pi Zi Z 2 Z 3 + Z 2 Za Z< + Za Z» Zi + Z» Zi Z 2 + Z j 

2pocWH 


_ (Zb + pc) ^ 2^ ^ 


Zb zj 


+ pc 


where 


zj=4+-^ 


ft A? 

i = 1 ^ 


(2.2.18) 

Section 2.2.1: Nondimensionalizatton 

Nondimensionalization of the general results presented 
above makes use of the nondimensional groupings shown in Table 
2.2.1 below. Typical numerical values are based on a square 
aluminum panel of length 12" on each side and 1/16" thick. Note 
that noncompensated mass ratios place the ART design frequencies 
(the frequencies of maximum noise reduction) at co = 1.0 (two 
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Table 2.2.1: Base nondimensional parameter components and values for the four panel and two 
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panel model) and a> = 1.0, 2.0, and 3.0 (four panel model) for the 
branch analysis without apparent mass. Compensated mass ratios 
place the ART design frequencies at the same frequencies but are 
compensated to include the effect of apparent mass loading, which 
will be discussed in detail later. Much of the analysis results 
presented in this work will use similar nondimensional groupings. 

Section 2.2.2; Branch Analysis Results 

Figure 2.2.4a presents the transmitted nondimensional 
pressure ratio pt/pi for both identical and ART tuned panels. For 
identical panels with a natural frequency of 1.0, note that the 
pressure ratio approaches a value of 1.0 at the natural frequency. 
(In the limit of no damping, this pressure ratio would be exactly 
1.0 at the natural frequency.) However, ART tuned panels show 
pressure ratios much less than 1.0 between the two panel natural 
frequencies of 0.5 and 1.5. Both results approach the same 
limiting values at very low and very high frequencies. Figure 
2.2.4b presents noise reduction results in decibels (dB) for the 
same configuration. Both identical and ART tuned panels show 
stiffness dominated noise reduction behavior (-6 dB/octave) at low 
frequencies and mass dominated noise reduction behavior (+6 
dB/octave) at high frequencies. 

As mentioned earlier, this technique provides simple results 
to check more complicated analyses. While apparent mass effects 
are not included (and are an important part of the general solution 
of the problem), the branch analysis represents the easiest 
analytical way to obtain quantitative results. 
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Nondimensional Frequency 

Figure 2.2.4a: Branch analysis pressure ratio results 

for identical panels and ART tuned panels. Values 
for the nondimensional parameters are as shown 
in Table 2.2.1. 



Nondimensional Frequency 


Figure 2.2.4b: Branch analysis noise reduction 

results for identical and ART tuned panels. 
Values for the nondimensional parameters are 
as shown in Table 2.2.1. 
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Section 2.3: The Full Four Panel Analysis 


This section presents the analysis of an idealized wall 
composed of multi-panel ART subsystems. The set of N panels in 
each subsystem grouping is assumed to repeat periodically along 
the paneled wall. As mentioned in Chapter 1, if each ART 
subsystem contains N panels tuned to N different resonance 
frequencies, there can occur N-l frequencies where the panels 
produce acoustic cutoff of the transmitted sound field. This fact 
demonstrates the possibility of canceling several propeller blade 
passage harmonics in an aircraft application, as demonstrated in 
Figures 1.2.1 and 1.2.2, for example. 

An infinite wall of panels composed of identical panel 
subsystems which repeat periodically in all directions is shown in 
Figure 2.1.1. The coordinate system is now located in the center of 
the panel grouping, and the panels are labeled in the matrix-like 
subscript notation as shown in Figure 2.1.1 and Figure 2.2.1. Note 
that this is the same analysis configuration used in the branch 
analysis described in Section 2.2.2, except for the coordinate 
system location. Suppose each panel subsystem contains four 
rectangular panels, with each panel having height H and width W, 
as shown in Figure 2.2.1. As before, each of the four panels can 
have different dynamic properties; however, other cases involving 
subsystems with one, two, or three panels can be recovered by 
defining various combinations of the panels to have identical 
dynamic properties. 

The analysis of normal incidence acoustic transmission 
through such a paneled wall has been accomplished using a 
number of analytical methods. The most efficient method first 
considers the forces associated with the motion of a single panel in 
each subsystem with the other panels held fixed. The 
corresponding forces associated with the other panels can be 
deduced from the symmetries of the problem. The general 
acoustic field for arbitrary panel motions can then be constructed 

26 



using superposition. The dynamics equation for each panel 
determines the relationship between the acoustic pressure on each 
panel and the resulting panel velocity. For simplicity, as in the 
branch analysis, the panel dynamics are modeled in terms of the 
impedance of an equivalent spring-mass-damper system as 
previously described by Equation (2.2.3). 

The acoustic pressure p(x,y,z)c * wt satisfies the reduced 
wave equation 


V 2 p + (co/c) 2 p = 0 (2.3.1) 

The acoustic velocity vector is [u(x,y,z)i + v(x,y,z)j+ w(x,y,z)k] e * wt . 
The x-component of the acoustic momentum equation is 

icou = - 1*— (2.3.2) 

P dx 

Therefore the velocity component u(x,y,z) can also be shown to 
satisfy the wave equation 

V 2 u + (co/c) 2 u = 0 (2.3.3) 

Now consider the motion of the panel denoted by “12” 
located in the first quadrant, as shown in Figure 2.2.1, with the 
other panels held fixed. The panel velocity mode shape is given by 
C>(y,z). With the other panels held fixed, the velocity boundary 
condition on the panel wall (x = 0) for the acoustic field is given by 

0( ) = |f(y* 2 ) forOSySWandOSz<Hl (2}4) 

1 0 for -W £ y £ 0 and/or -H £ z < 0 
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Since the panel subsystems are identical, <J>(y, z) = 0(y + 2W, z) - 
<D(y, z + 2H) = <D(y + 2W, z + 2H). The periodic function <J>(y, z) can 

be represented by a Fourier series of the form 

oo oo 

0(y, z) = aoo + S Za mn sinful sinM (2.3.5) 

For the present analysis, the panels are assumed to move 
as rigid spring-mass-damper panels; therefore, (j>(y,z) = constant = 
1.0*Ui2. The Fourier series for the velocity boundary condition 
then takes the form 


„ 12 (0,x, y ) = Ul2 + yi2£__l_sin^ii^ 

l2 ^ 07 4 * „.i(2n-l) W 


+ MUX_I sin (2m-l)* z 

* m=i (2m-l ) H 


(2.3.6) 


oo oo 


+ Hl2£ £ 


tc 2 m =i n=i (2n- 1 )(2m- 1 ) 


. (2n-l)jcy (2m-l)nz 

sm — — sin 

W H 


Since u(x,y,z) satisfies the wave equation subject to the boundary 
condition ui2(0,y,z) given above, the solution is 


Ul2(0 , x , y) = ^ + 


+ Mil £ _J sin 

* m=i (2m-l) H 


(2.3.7) 


oo oo 


—A sin ^"- 1 )** . sin < 2m - 1)ltz - 

tc^ m=i n=i (2n- 1 )(2m- 1 ) W H 


+ U12 £ £ 


The wave numbers in this expression are given by 
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(2.3.8a) 


koo=k = a 


kon = 


kmO ~ 



(2.3.8b) 


(2.3.8c) 


kmn — 


J 

(2m - 1) 7t 

2 

4_ 

(2n - 1) jc? fttT 

V 

H J 

T 

W J lcJ 


(2.3. 8d) 


Note that the panels are assumed to be acoustically compact and 
that only the one-dimensional mode propagates; thus ©H/c and 
toW/c < K. 

The x-component of the momentum equation can then be 
integrated to obtain the corresponding pressure solution 


(0,x,y) 


kopUi2 


= P2}hl e - iko»x + it»pUi2 T 1 sin (2n-l)»y . e -ko, 

4 ” n-i(2n-l)k 0 . W 

4M 

+ icopUi 2 ^ 1 sin ( 2 m-l)tc 2 g . km0 x 

* m=i (2m-l)k m o H (2.3.9) 

. X £ 4 S in sin (2m ' 1)7C --g' 

m=i n=i (2n-l)(2m-l )k mn W H 


The next step in the derivation is to determine the forces 
on each panel in the array due to the motion of panel 12 by 
integrating the pressure expression over each panel surface. Let 
Fjj kl be the force on panel ij due to the motion of panel kl. Then 


12 


< 


Pi 2 (0+>y.z)dyd z 


(2.3.10) 
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where A 12 indicates the surface of panel 12 (referring to Figure 
2.2.1); the factor of two in front of the integral accounts for 
radiation from both sides of the panel. Carrying out the 
integrations for the forces induced on each panel by the motion of 
panel 12 yields 

F 12 12 = 2pcWH[L + S n + S m + S mn 1 U 12 * C s U 12 (2.3.1 la) 

L4 J 

F n 12 = 2pcWH [-L - S n + S m - S mn lUi 2 = C w Ui 2 (2.3.11b) 

1.4 J 

F 22 12 = 2pcWH fi. + Sn - s m - S m „l U 12 = Ch U 12 (2.3. 1 lc) 

L4 J 

F 21 12 = 2pc WH \ L - Sn - s m + Sm„l U 12 = C D U 12 (2.3.1 Id) 

L4 J 


where 

00 

Sn = 2ik 2 1 (2.3.12a) 

7C 2 n=l (2n - l) 2 kon 


S m = 2ik X 1 (2.3.12b) 

7t 2 m =i (2m - l) 2 k m o 


= Mik X I 


2j L, 1 

Jt 4 n=im=i (2n - l) 2 (2m - l) 2 k mn 


Problem symmetry and superposition now allow these 
results to be generalized to the case where all panels are in motion. 
Let Fjj be the total force on panel ij due to the motion of all other 
panels. In general, it can be deduced that 
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(2.3.13) 


Fjj = Ui iC*i 2 -t- Ui2Cfj 2 + U 2 iC^j + U22C?22 


The coefficients CkiU are related as 


Cii = c ',2 = d 2 ',= c? 2 = Cs 

(2.3.14a) 

C?u= Ci 2 = c?,= $ 2 = Cw 

(2.3.14b) 

C?,= C?2= C* 2 2 . = C22= C D 

(2.3.14c) 

c?ii = Cj 2 = d 2 \= c 22 = c H 

(2.3. 14d) 

Therefore, the acoustic forces due to panel motion 
as 

can be written 

Fn = UnQs + U 12 CW + U 2 iC H + U 22 Cd 

(2.3.15a) 

Fi 2 = UnCw + Ui 2 C s + U 2 iCd + U 22 Ch 

(2.3.15b) 

F 2 i = UhCh + Ui 2 Cd + U2iCs + U 22 Cw 

(2.3.15c) 

F 22 = UhCd + Ui 2 Ch + U 2 iCw + U 22 CS 

(2.3. 16d) 


In the complete problem, acoustic forces also arise from a pressure 
wave striking the panel barrier at normal incidence. By the 
principle of superposition, the wall can be treated as rigid when 
including the effect of this incident pressure. The resulting wall 
motion, actually in response to this pressure wave, can be treated 
separately using the solution developed above. The appropriate 
wall response due to this forcing is then found by analyzing the 
dynamics of the panels; the latter is contained in the panel 
velocities. The right traveling incident wave, striking the panels at 
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x = 0- and reflecting as if from a rigid wall, has amplitude Pi. The 
corresponding pressure field is 2Pi cos(kx) e i®*, where the factor 
of 2 accounts for double amplitude due to the reflection from a 
perfectly rigid wall, and the resultant pressure is the real part of 
the complex pressure and harmonic in time. Note that adding the 
previous pressure field solution due to panel motion corrects the 
reflected pressure to its actual value and accounts for the 
transmitted wave. The force exerted by the incident wave and its 
perfect reflection on each panel is 2PiWH. The net acoustic force 
on a given panel is therefore 

F N ct ;: = Fij + 2PiWH (2.3.17) 

The above results are used in the force balance equations 
for the panel motion. Suppose each panel is modeled as a spring- 
mass-damper system of mass Mij connected to a rigid frame by a 
spring constant Sij and a damper Rjj. The net force and velocity of 
each panel are related through the mechanical impedance Zij as 

FNij = Zjj Uij = [ Rij + i (©My - ] Uij (2.3. 1 8) 

Writing out this equation for each of the four panels yields a 
system of coupled equations for the panel velocities 


(Zii +Cs) 

C W 

Ch 

Cb 

'Un" 


' 2 P 1 WH " 

C w 

(Zi 2 + Gs) 

Cb 

Ch 

Ui2 

, 

2PiWH 

Ch 

Cb 

(Z 2 i + C s ) 

Cw 

U 2 1 


2PiWH 

Cb 

Ch 

Cw 

(Z22 + C S ). 

-U 22 . 


. 2PiWH . 


(2.3.19) 
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In general, it is necessary to solve this equation numerically. 
However, one closed form result for the full problem will be 
mentioned later. 

Having solved the governing equations for the panel 
velocities, the pressures upstream and downstream can be 
computed and the noise reduction across the wall can be 
determined. It is readily shown by superimposing the earlier 
results for the motion of a single panel that the net transmitted 
one-dimensional pressure wave (x > 0) is given by 

Pl D (x,y,z) = P t e i (“ t ' kx > = — [Un + U12 + U 2 i + U22] e i(tDt ' kx) 

4 

(2.3.20) 

Similarly, the one-dimensional wave field on the incident side (x < 
0) is given by 

Pi. D (x,y,z) = 2Pi cos(kx) e l<ot - — [Un + U12 + U21 + U22] £ l(tot ‘ kx) 

4 


(2.3.21) 

These expressions can be used to determine the noise reduction 
across a paneled wall. Note that the above expressions do not 
include the higher pressure modes in the near field of the wall. 
These modes are less important from a noise standpoint because 
they decay exponentially and are nonpropagating; the effect of 
these modes has been included in the solution for the panel 
motions. 

As mentioned earlier, certain relatively simple closed form 
solutions can be obtained for special cases of the above analysis. 
One such case involves a panel subsystem with only two types of 
panels. Consider the particular case of two-dimensional panels of 
height H with unit width. After considerable work it is possible to 
show that the ratio of transmitted pressure to incident pressure 
takes the form 
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Pl 

Pi 


(Z e i + Ze2) 


(2.3.22) 


(Zd+Zrf) + ^ . if^Ze, + Ze 2 ) 
pcH 2pcH 


where Zei = Zi + io)m a and Ze 2 = Z 2 + i©m a are equivalent net 
mechanical impedances, with the impedances Zi and Z 2 being the 
mechanical impedances of the two types of panels. The quantity 
mj is analogous to an equivalent hydrodynamic mass associated 
with the higher order evanescent pressure modes produced by the 
out-of-phase components of panel motion, and may be expressed 
as 


m a = 


_L 


16pcH 


7t 2 co n = i(2n - 1) 


, [2n-l]rcc/ 

2o>H 


1- V 2 


- 1 


(2.3.23) 


Note that Equation (2.3.22), with m a = 0 recovers the branch 
analysis result for a two panel system. The four panel model will 
effectively yield this result when neighboring panels are designed 
with the same dynamical properties. 
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Section 2.4; Analytical Results 


Section 2.4.1; Two Panel R fiSMltt and Discussion 

Figure 2.4.1 shows a noise reduction prediction in decibels 
plotted against nondimensional frequency for identical and ART 
tuned panels. Parameters are as shown in Table 2.2.1 under 
“Noncompensated 2 Panel Model”. (Panel mass compensation is 
addressed later in this section.) These parameters correspond to 
square aluminum panels measuring 12" on a side, having a 
thickness of 1/16". As with the branch analysis, both panel 
systems show similar stiffness dominated noise reduction behavior 
at low frequencies (-6 dB/octave) and mass dominated noise 
reduction behavior (+6 dB/octave) at high frequencies. At a panel 
natural frequency, that particular panel will appear acoustically 
transparent in the limit of zero damping. Between the frequencies 
of 0.8 and 1.2, the ART contribution to noise reduction is greater 
than 20 dB. This result differs only slightly from the branch 
analysis result in that the peak ART noise reduction frequency has 
decreased slightly, as shown in more detail in Figure 2.4.2. This is 
due to the additional mass loading of the fluid surrounding the 
panels. From a practical viewpoint, the frequency regime around co 
= 1.0 ± 0.2 indicates a fairly large noise reduction bandwidth; that 
is, a frequency regime where the noise reduction is greater than a 
certain amount (in this example, about 20 dB). This fact may be of 
value in reducing large sound levels at slightly varying 
frequencies; for example, a slightly varying propeller blade 
passage frequency due to propeller gust loading. The large noise 
reduction bandwidth phenomenon also indicates the robust nature 
of this tuning method. 

The solution of the governing equations for panel complex 
velocities permits the extraction of phase data for the individual 
panels. Figure 2.4.3 shows the theoretical phase difference 
between the two ART tuned panels including apparent mass as a 
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Nondimensional Frequency 

Figure 2.4.1: ART and identical two panel noise 

reduction prediction with apparent mass plotted 
versus nondimensional frequency. Values for 
the nondimensional parameters are as shown in 
Table 2.2.1. 



Nondimensional Frequency 

Figure 2.4.2: Detail of the branch analysis and 

full theory noise reduction predictions around 
the ART cancellation frequency. Remainder of 
noise reduction calculation is essentially the same. 
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function of nondimensional frequency. Here, the damping ratio 
has been fixed to 0.01, and as a result of this low damping, much 
of the frequency domain between the natural frequencies 0.5 < co < 
1.5 shows almost perfect out-of -phase behavior. Figure 2.4.4 
shows panel velocity magnitude relationships, where the velocity 
magnitude is relative to a constant amplitude of incident pressure 
Pi. (This differs qualitatively from the concept of Figure 1.3.1, 
where panel velocities are sketched relative to constant forcing.) 
The panel velocity magnitudes are equivalent at just below co = 1.0, 
which corresponds to the “ART design frequency”, or co ART - 
However, as shown by Figure 2.4.1, significant noise reduction can 
occur at nondimensional frequencies where the panel velocity 
magnitudes are somewhat different. Figure 2.4.4 indicates that the 
panel motions are coupled to the surrounding pressure field in a 
frequency dependent manner. Note that at the individual panel 
resonances, the non-resonance panel velocity is at a minimum. 
Again, at a natural frequency, a panel will appear almost 
acoustically transparent, and nearly all of the sound transmission 
will occur through that panel. 

Section 2,4.2; Averaged Noninteracting Noise — Reduct i on 

In order to observe the true contribution of out-of-phase 
motion to noise reduction across the panel configuration, it is 
necessary to compare the ART data to an idealized case where only 
uniform panels are used, and the resulting acoustic fields do not 
interact, as shown in Figure 2.4.5. This comparison allows the 
mass law and stiffness law effects associated with moving the 
panel resonance frequencies apart to be separated from the effects 
of acoustic cancellation or cutoff which are central to the ART 
concept. Such an averaged noninteracting noise reduction (ANNR) 
may be derived for N panels as 

(NR) A nnr(“) = 10 log(N) - £ 1 0' (NR,i< “ ) /i 0 (2.4.1) 

i = 1 
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Nondimensional Frequency 

Figure 2.4.3: Theoretical phase difference in 

degrees between the ART panels with apparent 
mass as a function of nondimensional frequency. 



Figure 2.4.4: Panel velocity relationships from 

two panel theory with apparent mass. 
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Figure 2.4.5: Averaged Noninteracting 

Noise Reduction (ANNR) schematic. 


Figure 2.4.6 shows the two panel ART noise reduction plotted with 
the ANNR result. Note that the ANNR result is intuitively 
satisfying, since its behavior is stiffness and mass dominated at 
low and high frequencies. ANNR also reaches a minimum at 
resonance frequencies, and this noise reduction minimum for 
identical panels is slightly better than the ART minimum since 
ANNR has another set of panels (the higher resonance frequency 
set) contributing to the noise reduction effort. However, ART still 
indicates a true maximum noise reduction increase of 20 dB or 
more between the panel natural frequencies, depending on the 
desired size of the noise reduction bandwidth. 

Section 2.4.3: Four Panel Results and Discussion 

Figure 2.4.7a shows the pressure ratio pt/pi for both 
identical panels tuned to ton = 1.0 and ART panels tuned to co n = 
0.5, 1.5, 2.5, and 3.5. Other parameters are as listed in Table 2.2.1. 
All previous comments with respect to pressure maxima and 
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Figure 2.4.6: Two panel ART noise reduction plotted with 

Averaged Noninteracting Noise Reduction (ANNR) with 
apparent mass. 


minima, mass and stiffness domination, and acoustic transparency 
apply here. However, Figure 2.4.7a demonstrates the ability of an 
ART tuned panel array with four different panel tunings to cancel 
three frequencies. (Measurements by Woodward on a scale model 
SR-7A high speed propeller have shown that the tonal content of 
the propeller noise spectra was typically limited to the first three 
harmonics, with higher tone orders either not present or masked 
by broadband background noise probably produced by the flow 
over the microphone.^) Figure 2.4.7b shows the equivalent data 
as a noise reduction in decibels. Again, note the fairly wide 
bandwidth of noise reduction at 20 dB, indicating the robustness of 
the tuning method. With four panels, the apparent mass loading 
becomes more important, especially with respect to the higher 
natural frequency panel since this panel has the lowest mass. The 
lowest ART cancellation frequency displays the greatest noise 
reduction simply because these lower frequencies are more cutoff 
due to their longer wavelengths; restated, the panels are more 
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Nondimensiona! Frequency 

Figure 2.4.7 a: Four panel ART pressure ratio P t / Pj 

plotted against identical panels pressure ratio with 
apparent mass. 



Nondimensional Frequency 


Figure 2.4.7b: Four panel ART and identical panels 

noise reduction predictions with apparent mass. 
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acoustically compact at lower frequencies. Figure 2.4.8 shows the 
same ART configuration plotted with the ANNR results for four 
panels, showing a true maximum gain in noise reduction of at least 
15 dB at io ART = 3.0 and at least 20 dB at ©art = 10 over the 
identically tuned system. 

Section 2.4.4; Mass Ratio Optimization 

Figure 2.4.9 shows the ART four panel model results plotted with 
and without the apparent mass loadings. Note that the data 
without apparent mass is equivalent to the branch analysis results 
discussed earlier and presented in closed form via Equation 
(2.2.18). Both calculations have input data which should set the 
ART cancellation frequencies at co A RT = 1A 2.0, and 3.0. However, 
in the case of the simulation with apparent mass, the additional 
loading of the fluid surrounding the panel adds mass, and thereby 
reduces the natural frequencies of the combined panel/fluid 
system. It is generally desired to place the maximum ART 
cancellation frequencies at specific values in order to attenuate 
harmonics of the fundamental frequency, or to be able to control 
precisely the frequencies at which noise reduction will occur. This 
can be accomplished by altering the panel natural frequency, 
stiffness, or mass. Mass ratio (that is, the ratio of any panel mass 
to a reference panel mass in the panel array) optimization is 
chosen because this method is successful in laboratory simulations. 
Iteration of the four panel theory via a multi-dimensional Newton- 
Raphson method with equations of the form 

(o) ART - CD ART ) = 

V desired current/ 


Aa) ARTi 



Ei 

£i 


l d(MR)j 


d(MR)j 


(2.4.2) 
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Nondimcnsional Frequency 


Figure 2.4.8: Four panel ART noise reduction 

prediction with apparent mass plotted against 
Averaged Noninteracting Noise Reduction (ANNR). 



Nondimensional Frequency 


Figure 2.4.9: Noise reduction prediction using the ART 

four panel theory with and without apparent mass 
modes as indicated above. 
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yields three equations with three unknowns, one for each mass 
ratio. The mass of panel 11 is held fixed at 1.0, and the numerical 
solutions to Equation (2.4.2) are the adjustments to the current 
mass ratios MR,j to bring the ART maximum cancellation 
frequencies closer to the desired value. Equation (2.4.2) is iterated 
until convergence within 0.001 nondimensional frequency is met. 

Of particular interest in this problem is the efficient 
location of a minimum of a function within specified bracketing 
limits. The four panel pressure ratio plotted in Figure 2.4.7a shows 
an approximately parabolic and well-behaved function which is at 
a minimum when noise reduction is at a maximum. Brent s 
method provides an effective and quick scheme for finding minima 
in such situations. 19 Parabolic interpolation using three points is 
attempted first, using the panel natural frequencies which bound 
the minimum and a midpoint between them as initial guesses. If 
the parabolic steps become nonconvergent in finding a minimum, 
the routine switches to a golden section search. In the worst case, 
the routine alternates between these two methods until a suitable 
minimum is found. In practice, this method found a pressure ratio 
minimum between panel natural frequencies in about five to ten 
attempts, and generally found the optimum mass ratios after three 
to five iterations of Equation (2.4.2). The initial guesses for the 
mass ratio values were obtained using MACSYMA to solve three of 
the four Equations (2.3.19) with Cs = Cw = Oh = Cd - 0.0, an< J 
freezing the mass ratio of panel 11 at 1.0. 

Figure 2.4.10a shows the full four panel theory with and 
without compensated mass ratios. Note that the peak ART noise 
reduction values align to coart = 1*0, 2.0, and 3.0, and the peak 
noise reduction values increase. To further show the sensitivity of 
the four panel theory on panel size and apparent mass, Figure 
2.4.10b shows noise reduction prediction data for panels similar to 
the type used in the laboratory confirmation experiments (see 
Chapter 5). Here, the panel is assumed to be square and 4" on each 
side, and the panel mass is approximately 1 gram. Apparent mass 
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Nondimensional Frequency 

Figure 2.4.10a: Four panel theory with and without 

compensated mass ratios. Note shift of peak ART 
cancellation frequencies and slight increase in noise 
reduction for compensated system. Values are as shown 
in Table 2.2.1. 



Nondimensional Frequency 


Figure 2.4.10b: Full four panel theory with parameters 

changed as shown above. Other parameters equivalent 
to Figure 2.4.10a. 
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effects are consequently much more evident, as the 
uncompensated ART data shows that the peak noise reduction 
frequencies have severely dropped due to a higher apparent 
mass/panel mass ratio. Figure 2.4.10c shows a comparison of the 
values of the mass ratios for the previous two figures. Both 
noncompensated noise predictions begin with the same panel mass 
ratios. The compensated mass ratio values are then compared in 
the legend. As the panels with the higher natural frequencies and 
smaller size are less massive, the mass ratio optimization 
procedure compensates by adding more mass to the smaller 
and/or lighter panels. 

Finally, Figure 2.4.11 shows the magnitude of the 
nondimensionalized apparent mass sums as a function of 
nondimensional frequency for the standard parameters as shown 
in Table 2.2.1. Note that panels horizontally adjacent and 
vertically adjacent have the same apparent mass contributions to 
the governing Equations (2.3.19); that is, Cw = Qt- Equivalently, 
from Equation (2.3.12), S n and S m are identical. This again is a 
consequence of the symmetry of the problem. The apparent mass 
sums shown in Figure 2.4.11 are the coupling terms appearing on 
the off-diagonal elements of the linear system denoted by 
Equation (2.3.19); their contribution to the full problem by 
coupling panel motions to the surrounding pressure field increases 
as frequency increases. Also, note that adjacent panels produce 
stronger coupling through the apparent mass than diagonal panels. 
This phenomenon is readily explained by the perimeters matching 
along a line for adjacent panel boundaries and only matching at a 
point for diagonal panels. Adjacent panels slosh fluid from each to 
the other, and are thus more highly loaded than panels meeting at 
just a point. In the latter case, much less surrounding fluid is 
shared by diagonal panels in the sloshing motion. 
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Mass Ratio Values 


■ all panels 

noncompensated 
H large panel 
compensated 
E3 small panel 
compensated 



to =0.5 to =1.5 w=2. 5 o) =3.5 

n n n n 


Figure 2.4.10c: Comparison of mass ratios for noncompensated 

and compensated predictions with varying panel size and 
apparent mass loading. 
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1 


panels diagonal 



Figure 2.4.11: Magnitude of apparent mass sums 

versus nondimensional frequency showing increased 
coupling at higher frequencies. 


Section 2.4.5i Model Parameter Studies 

Having established the basic acoustic behavior of the four 
panel system with identical and ART tuned panels, parameter 
studies with the ART tuned system can establish performance 
expectations over a wide range of conditions. Figure 2.4.12 shows 
the effect of varying panel damping on the ART model. Increasing 
panel damping ratios reduces the peak ART noise reduction values 
and increases the minimum noise reduction at the panel 
resonances. In effect, adding more damping makes the system less 
acoustically transparent at resonances. However, Figure 2.4.12 
clearly shows that a comfortable noise reduction bandwidth of 
about 20 dB can still be maintained around the ART design 
frequencies ±0.2 nondimensional frequency. This is again an 
indication of the robust nature of the method. 
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Nondimensional Frequency 

Figure 2.4.12: Parameter studies with the four panel 

model varying the panel damping ratios as indicated 
above. Nonadjusted parameter values remain as 
listed in Table 2.2.1. 



Nondimensional Frequency 


Figure 2.4.13: Parameter studies with the four panel 

model varying sound propagation speed via the 
parameter co r H/c. The effective sound speeds are 

shown in the legend. 
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Figure 2.4.13 shows the effect of the variation of sound 
speed on the model. Increasing sound speed tends to reduce cutoff 
behavior and causes the ART cancellation peaks to become 
sharper. It is similar to variation of panel damping also in that 
minimum noise reductions at panel resonances decrease. Note that 
this is only an exercise in variation of a parameter in the computer 
code; physical panel systems under atmospheric loadings would 
not experience this degree of sound speed variation. It does 
indicate that different behavior might be expected if this tuning 
scheme were to be used in water, for example. However, the latter 
scenario would drastically increase the fluid loadings on the 
panels, and would require different analyses. 

Figure 2.4.14a shows variation of the panel aspect ratio 
while holding apparent mass loading fixed. This is clearly a higher 
order effect, as shown in Figure 2.4.14b, which shows a detail of 
the frequency regime around the highest ART frequency. This 
behavior shows that as the panel becomes smaller, apparent mass 
loading plays only a small effect on the higher natural frequency 

panel, as shown in Figure 2.4.14b. 

Figures 2.4.15a and 2.4.15b show more of the coupling 
between apparent mass and panel aspect ratio. Figure 2.4.15a 
shows variations in apparent mass loading by varying the density 
of the surrounding fluid from normal air density to eight times 
normal air density with a panel aspect ratio of 1.0. As can be 
expected, the less massive panels are more affected by the 
increased loading, and the noise reduction bandwidths begin to 
narrow with higher fluid loadings. The prediction for eight times 
normal air loading in Figure 2.4.15a even shows the emergence of 
an additional noise reduction peak at a frequency slightly above 
the highest ART cancellation peak. This might be viewed as a 
"spring-like" inertance effect of the fluid around the panel. Figure 
2.4.15b shows the same variations in apparent mass loading on 
panels with W/H = 0.25, and the inertance effect is less prominent. 
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Nondimensional Frequency 

Figure 2.4.14a: Parameter studies with the four 

panel model varying the panel aspect ratio W/H 
as shown above. Effective panel widths are shown 
in the legend. Panel height remains at 12”. 



Figure 2.4.14b: Detail of highest ART cancellation 

frequency increase due to panel aspect ratio change. 
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Nondimensional Frequency 

Figure 2.4.15a: Parameter studies with the four 

panel model varying the apparent mass loading 
parameter as indicated above. Panel aspect ratio 
W/H = 1.0. 



Nondimensional Frequency 

Figure 2.4.15b: Parameter studies with the four 

panel model varying the apparent mass loading 
parameter as indicated above. Panel aspect ratio 
W/H = 0.25. 
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Figure 2.4.16: Parameter studies with the two panel 

model varying the higher panel natural frequency 
as indicated above to increase the noise reduction 
bandwidth. Nonadjusted parameter values remain 
as shown in Table 2.2.1. 


Figure 2.4.16 shows how spreading the panel natural 
frequencies apart can effectively increase the noise reduction 
bandwidth. This application may hold promise for noise problems 
which are both narrowband and broadband, as mentioned by 
Knowles. 20 

Figures 2.4.17a and b show the effect of variation in the 
termination impedance in the four panel model. Figure 2.4.17a 
varies the real portion of the termination impedance from anechoic 
(Zb = 1.0 + iO) to (Zb = 8.0 + iO) by successive doubling. Note that 
the duct is 10 panel widths long and the microphone is placed at a 
position 7 panel widths downstream. The noise reduction 
predictions now become contaminated by the reflected wave. 
However, the ART noise reduction peaks clearly remain intact, 
although the effect of “hardbox” frequencies can be seen, 
intertwined with the minima associated with panel natural 
frequencies. These hardbox frequencies are denoted by 21 
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(2.4.3) 


f m = 


me 

2Lduct 


where m is an integer. In nondimensional terms. Equation (2.4.3) 
becomes 


% 


mc- 

S1S4 


(2.4.4) 


For the parameters shown in Table 2.2.1, this means that noise 
reduction minima will occur at frequency multiples of 0.561 (as 
well as at panel natural frequencies), and close examination of 
Figure 2.4.17a shows this to be the case. However, an interesting 
feature of the ART method is its ability to work at a hardbox 
frequency; examination of both Figure 2.4.17a and b shows that 
while the hardbox frequency at co = 2.80 is visible, is still 
effectively attenuated in the surrounding frequency regime by the 
ART concept. Additionally, Figure 2.4.17b shows wider variations 
in the downstream termination impedance, from (Zb = 1.0 + iO) to 
(Zb = 50.0 + iO). However, any termination impedance equal to or 
harder than (Zb = 10.0 + iO) looks essentially the same. However, 
care should still be exercised in this regard; placing a hardbox 
frequency exactly at an ART design frequency in a single panel 
wall-duct configuration will deteriorate performance 
unnecessarily. 

Noise attenuation has been shown to be enhanced by 
increasing the panel damping in a structure.22 Figure 2.4.18a 
shows that while the ART concept is more effective at lower values 
of the damping constant £, it can still contribute to noise reduction 
in structures with higher damping. The noise reduction around the 
panel natural frequency regime ± 0.2o> has remained nearly the 
same despite increased damping. Figure 2.4.18b compares a very 
highly damped ART noise prediction (C = 0.4) to the ANNR noise 
prediction for uniform panels with the same damping; ART still 
outperforms uniform panels, as can clearly be seen. Note that the 
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Nondimensional Frequency 

Figure 2.4.17a: Parameter studies with the four panel 

model varying the real part of the nondimensional 
termination impedance as indicated above. 



Nondimensional Frequency 

Figure 2.4.17b: Parameter studies with the four panel 

model varying the real part of the nondimensional 
termination impedance 7 ^ as indicated above. 
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Nondimensional Frequency 


Figure 2.4.18a: Parameter studies with the four panel 

model for high values of damping as indicated. 



Nondimensional Frequency 

Figure 2.4.18b: High damping comparison (C=0.4) 
of ANNR and ART noise reduction predictions. 
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lowest frequency portion of the ANNR calculation occurs at the 
lowest panel natural frequency co = 0.5, another indication that low 
frequencies are the most difficult to attenuate. Figure 2.4.18b 
suggests that the ART concept might also be integrated with other 
noise reduction techniques such as increased damping. 

Finally, Figures 2.4.19a and 2.4.19b show a comparison 
between tuning configurations as indicated in Figure 2.1.2. The 
data was generated by moving panel natural frequencies with 
their associated mass ratios to match the possible tuning patterns 
from Figure 2.1.2. Note that when apparent mass is compensated 
using panel mass ratio adjustment, the different systems behave in 
remarkably the same way. Only the highest ART cancellation 
frequency (due to the combination of lowest mass panel with the 
highest natural frequency) is affected in a noticeable way. Minor 
mass ratio compensation could clearly be done on any 
configuration, and effectively locate the noise reduction 
attenuation maximums at any frequency region desired. 
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Nondimensional Frequency 


Figure 2.4.19a: Comparison of tuning methodologies 

for the four panel model using the 3 patterns 
indicated in Figure 2.1.2. 



Figure 2.4.19b: Enlargement of noise reduction 

results around ART cancellation nondimensional 
frequency = 3.0 to highlight subtle differences 
between tuning strategies. Patterns as in 
Figure 2.1.2. 5 8 



Chanter 3 

External Pre ssure Field Modelin£ 


Section 3.1; Introduction 

The purpose of this investigation is to determine the 
effectiveness of the ART concept under an external propagating 
pressure field such as that which might be associated with 
propeller passage by an aircraft fuselage. Figure 3.1.1 shows a 
schematic representation of a typical aircraft cross section near the 
propeller plane. Associated with the motion of the propeller near 
the fuselage is a pressure wave which sweeps by the fuselage 
panels at some frequency and wavelength which is related to the 
engine operating speed and propeller characteristics. (Nallasamy, 
et. al., have shown experimentally that the measured acoustic 
pressure signal in the plane of the propeller is very closely 
sinusoidal at a number of blade pitch angles. 23 ) In general, the 
problem deals with modeling the interaction of fluid and structural 
components within the realm of ART tuning. 


Section 3.2: Analysis Derivation 

The physical schematic for the system is presented in Figure 
3.2.1. The analysis configuration consists of a two-dimensional 
duct, where the position of the coordinate system is shown in the 
lower left hand side of the schematic. A linear array of Q panels, 
each of height H, is placed along the z axis at the location x = 0. For 
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Figure 3.1.1: Front view of aircraft fuselage showing 
relative locations of propeller and panels in the propeller 
plane. Dimensions for aircraft used in computational cases 
are given in Table 3. 1 . 1 . 
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Direction of x 
External Pressure 
Field Motion 
(positive upward) 

Figure 3.2.1: External pressure field modeling schematic 
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convenience, these panels are modeled as standard acoustic 
spring-mass-damper impedances. Note that these Q panels may be 
tuned in any manner. For example, a standard ART tuning of 2 
panels per ART panel pair may be used (one panel has a low 
natural frequency, and the second panel has a high natural 
frequency) or more complex tuning arrangements may be used 
where n panels have n different tunings, where the only constraint 
is n < Q. The choice for the panel tuning methodology is not 
required at this point in the analysis. The termination impedance 
for the duct is considered to be anechoic. A propagating external 
pressure is present on the exterior side of the panel wall barrier 
(at x = 0-) and is of the form 

Pex = P e <? i(<Bt -k z) (3.2.1) 

This pressure is assumed to be sweeping past the panel array, as 
shown in Figure 3.2.1. 

The analysis begins by considering that one panel of the Q 
total panels is raised on both the duct and an associated image 
duct. For the duct. 


f(z) = 0, QH > z ;> rH 

f(z) = 1 , rH > z > (r-l)H (3.2.2) 

f(z) = 0, (r-l)H > z £ 0 

where r is a panel counting integer 1 i r < Q. Similarly, for the 
image duct 

f(z) = 0, 0 £ z 2: -(r-l)H 
f(z) = 1, -(r-l)H £ z£ -rH (3.2.3) 

f(z) = 0, -rH > z 2: -QH 

A Fourier series representation of these raised panels on both the 
duct and the image is then given by 
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fW = Q + n | C ° S W ' 2)) Si " fell C0$ ftp) (3 2 ' 4) 

From hence forth, one need not consider the image duct in the 
analysis. The velocity representation for the raised panel on the 
duct is then assumed to be 

u(0,z) = 0, QH £ z £ rH 

u(0,z) = U r e io>t , rH £ z S> (r-l)H (3.2.5) 

u(0,z) = 0, (r-l)H ;> z ;> 0 


Using the Fourier series representation in Equation (3.2.4) above, 
the general velocity Ur for panel r is then compactly written as 


u r (0,z,t) = U T e'°> 1 



( 3 . 2 . 6 ) 


The modal velocities u n (x,z) in the duct satisfy the linearized two- 
dimensional wave equation 


§3!. + ?!H2. + k 2 Uo =0 (3.2.7) 

dx 2 dz 2 

where 

uq = giko* — > ko (3.2.8) 


and 



(3.2.9) 
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Note that the mode zero velocity, denoted by Equation (3.2.8), 
always propagates. The higher modal velocities, denoted by Un, 
contain the wavenumbers k n ; some of the k n are cut off, while the 
remainder propagate. The frequency of cut on for the n 1 * 1 mode 
may be described as the solution to 


G)=c/U£- 


( 3 . 2 . 10 ) 


Similarly, at any given to, the maximum propagating mode will be 
denoted by the numeric truncation of 


coQH 
n = - - 

itc 


(3.2.11) 


Mode numbers up to and including this truncated value for n will 
propagate. All higher modes will be cut-off. For convenience in 
the remainder of the derivation, the truncated value of n 
represented by Equation (3.2.11) will be denoted as nlow. 

Similarly, an integer value of nhigh will be defined where nhigh = 
nlow + 1. 

The velocity contribution from the r th single moving panel at 
any point in the duct is then denoted by 


+IV 


U r (x,Z,t) = U r e ia>t ^e* ikoX 


nlow I 

i<ot | V e -ik B x j_4_ cos 

" InTt 

n=l \ 

I MODES-1 


W r - 2)) sin (£)) cos P 


+U ' Ci “‘ nl, h e l^ 0S W r ' 2)^" (“r P 


(3.2.12) 
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In this formulation, the modal wavenumber k n is forced to be a 
positive real number at all times. Therefore, for propagating 
modes, 



(3.2.13) 


while for cutoff modes, 


k c 

K n 


/nil 2 - 


\CH 


(3.2.14) 


Application of the momentum equation in the x direction 

p(x,z,t) = -itop | u(x,z,t)dx (3.2.15) 

results in an expression denoting the pressure contribution by the 
r th moving panel at any location in the duct. Specifically, 

p r (x,z,t) = 

Qc 


icot| n V W ^P |_4_ cos 

nrc 


+ V,e‘“' | I -je 

, n=l k n 


(‘sKi-esHS) 


nr fa - y ... ™ (“) 

\ n=nhigh K n ' 


(3.2.16) 


Only the cut on pressure modes will contribute to the overall 
complex acoustic pressure far down the duct. In the panel array 
near-field, however, the resultant pressure is a combination of the 
cutoff evanescent modes and propagating modes. 
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Next, an expression is required to obtain the forces on each 
panel in the panel array. Conceptually, panel r will cause a force 
on panel s, where s is another panel counting integer operating 
over the same limits as r. The general sign convention shall be 
that a positive panel velocity will effect a negative panel force. 
Additionally, the effect of the external sweeping pressure must 
also be included. This may be expressed as 

F rs = -A| p(x=0+,z,t)dz + I Pexdz (3.2.17) 

Jpu*\ limits 

The variable A = 1 to account for hydrodynamic loading on 
only one side of the panel array; A = 2 allows for hydrodynamic 
contributions on both sides of the panels. Setting A = 2 and 
allowing only one acoustic mode to be present (the one 
dimensional mode) along with a long wavelength for the external 
pressure forcing allows the recovery of the branch analysis results 
presented in Chapter 2 as a check of both analysis and coding 
accuracy. If the external sweep forcing is assumed much greater 
than the hydrodynamic loading, it can be expected that the 
solutions for A = 1 or 2 will not differ greatly except at the panel 
natural frequency (in this general frequency regime the panel 
velocities are the greatest) and perhaps at the ART frequency, 
where the out-of-phase motion of the ART panel pairs causes a 
different hydrodynamic loading on the panels. Allowing x = 0+ in 
Equation (3.2.16) above and integrating along individual panel 
limits defined by s yields 
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(3.2.18) 


where use has been made of a standard trigonometric identity to 
simplify the expression. Note that the integration described above 
is a piecewise application of orthogonality. Next, the integration of 
the pressure term shown in Equation (3.2.1) is also considered 
over each panel in the array as 

J r rH 

p e gi(©t - k«z) _ (g-ik.rH (j . (3.2.19) 

(r-l)H ke 


The following definitions may now be made for notational 
convenience. 


B 


prop 


nlow _ 

Y 1 - 

n=l k„ P 


JL 


cos 


\(n *f 




S - 

2/ 


sin 


BJL 

m 


( 3 . 2 . 20 ) 


B 


cut 


MODES- 1 „ 

Y -L 

" , C 

n=nhigh K 


i cos (^l r ^)) cas (^^)) 
' ' 


sm 


\ 

I 


(3.2.21) 
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The governing equations may then be expressed as 

iPt (^-ik.rH (] . c ikjl)) = 

k e 

Ur(A(pcH + copQH 2 Bprop + icopQH 2 Be*)) (3.2.22) 

+ U T Zr (iff r = s) 

where the common factor of e * <ot has been cancelled. As stated 
above in Equation (3.2.22), the term with velocity U r multiplied by 
the panel mechanical impedance term Z r is only added into the 
equation if r = s; i.e., if this coefficient in the matrix is on a diagonal 
term. Appropriate nondimensionalization then leads to a set of Q 
linear algebraic equations in the Q unknowns Ur(to). When these 
U r (co) are known, they may be reinserted into Equation (3.2.16) 
above (also in the appropriately nondimensionalized form) to 
obtain the pressure at any point in the duct as a result of the 
motion of panel r. The total overall acoustic pressure at any 
location in the duct is then obtained simply by summing the Q 
contributions of Equation (3.2.16). 

Section 3.2.1: Nbise Reduction Calculation 

Pressures on each side of the panel barrier contribute to the 
calculation of a nondimensional pressure ratio for the expression of 
a noise reduction in decibels at any location in the duct. Upstream 
incident pressure is a combination of the external sinusoidal 
forcing and the resultant loading from the apparent mass modes. 

Pa = Pex( z »0 + Pr( x - 0-. z »0 (3.2.23) 

= Pex( z »0 * Pr( x = 0+,z,t) 


Transmitted pressure pu is obtained by summing the Q 
contributions as given in Equation (3.2.16); together, pa and pb 
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form a pressure ratio of incident to transmitted pressure. The 
noise reduction in decibels is then denoted by 


NR d B = -20 log 


£b_ 

Pa 


(3.2.24) 


Section 3.2.2: Nftndimensiftnalizati<m 

Nondimensionalization of the general results presented 
above makes use of the nondimensional groupings as shown in 
Table 2.2.1 in Chapter 2. However, for this problem, a 
nondimensional velocity is defined as 

— _ UrppC (3.2.25) 

' Pe 

and a ratio of panel height to external wavelength nondimensional 
variable is defined as 


S 6 = K e H=22LH (3.2.26) 

As in Chapter 2, the standard set of parameters considers a square 
aluminum panel 12" on each side. More discussion on the code 
input parameters is given in section 3.3.1. 

Section 3.2.3; Numerical Solution Procedmg 

The numerical solution for the external pressure field 
modeling problem is accomplished in a two step procedure. The 
governing Equation (3.2.22) and its components are programed to 
allow for any number of panels to be included in the duct. 
Additionally, any number of panel impedances may be modeled in 
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any order, the latter being made possible by a definition of panel 
properties in a data input file. For a given value of 
nondimensional frequency to, the arrays B pr op an d B C ut are 
calculated. (Note that because these arrays contains the 
propagating and cutoff wavenumbers. Equations (3.2.20) and 
(3.2.21) are functions of co.) Equation (3.2.22) is first solved for the 
complex nondimensional velocities Ur, and these velocities are 
written to a data file. A second program then calculates pressures 
according to Equation (3.2.16), and the total nondimensional 
pressure, as mentioned earlier, is the sum total of all contributions 
of Equation (3.2.16). This procedure was adopted to allow multiple 
calculations of nondimensional pressures and noise reductions in 
decibels at many locations in the duct without recalculating 
velocities; the velocity calculation generally is the most time 
consuming portion of the procedure. Knowledge of the panel 
velocities for a given set of input parameters (panel size, mass, 
natural frequencies, apparent mass loading, damping ratios, and 
external propagating pressure speed) allows the calculation of 
pressures and noise reduction levels at any location in the duct. 

Two other programming issues are deserving of a short note. 
The numeric solution of the governing Equations (3.2.22) is first 
solved using LINPACK routines, which are available in the public 
domain; these routines are quite robust and check for matrix ill- 
conditioning. Given stable solutions for a number of cases, a faster 
matrix inversion routine was used; since the results obtained are 
identical, the faster solution routine is used for production runs. 
Also, given the desire to generate large amounts of data for a 
variety of panel configurations, a simple subroutine keeps track of 
all data inputs and file names for a given simulation. In this 
manner, all numeric “housekeeping” was computerized and no 
longer prone to human error. 


70 



Section 3.2.4i Results Verification 


The computer models mentioned in Chapter 2 may be used 
as a base upon which comparisons for ART and identical panels 
may be made. In particular, solving the sweep problem 
numerically for two ART panels with 1 acoustic mode and KeH -> 0 
will recover the ART branch analysis results. (Note that KeH 
cannot be set to exactly zero in the computer code; due to the 
formulation of the governing equations for the program, division 
by zero will result.) Additionally, the external pressure modeling 
computer program with many acoustic modes as input will be 
similar to the full 2 panel ART solution with apparent mass, except 
for results near frequencies where mode turn-on occurs. Also, this 
latter phenomenon can still be checked with an identical panels 
simulation using both the full 2 or 4 panel codes, as identical 
panels in a duct will produce the same noise reduction and 
nondimensional pressure results with and without apparent mass. 
Verification of the sweep problem results in this manner assures 
correct solutions and generates user confidence in the analysis and 
associated programs. 



Section 3.3; Analytical Results 


Section 3.3.1; Oase Study Inpuls 


Data from two aircraft are used in case studies for the 
sweeping analysis. Both simulations include the use of eight 
aluminum panels 12" on each side to construct the total height for 
the duct (fuselage sidewall) model. Table 3.3.1 shows parameters 
representing the two aircraft; the first is a Gulfstream Aerospace 
Commander 695 A aircraft, a high-wing business aircraft, as 
reported by Wilby and Wilby. 24 These input parameters will be 
referred to as the “long wavelength” case, due principally to the 
use of a three-bladed propeller. Additionally, this three-bladed 
propeller creates a low blade passage frequency, which is 
represented in the nondimensional parameter Si as shown in Table 
3.3.1. Note here that the external sweep wavelength is about two 
and one-half times longer than the overall duct height. Table 3.3.1 
also displays data from Kuntz and Prydz 25 representing a 
Gulfstream II aircraft modified for noise studies; this configuration 
will be referred to as the “short wavelength” case due to use of an 
eight-bladed propeller. The blade passage frequency of the latter 
case is correspondingly higher, as shown by the value of Si in 
Table 3.3.1 under the “short wavelength” category. In the latter 
example, the external wavelength is comparable to the duct height. 
The adjustment of the nondimensional parameter Si as shown 
places the lowest desired ART cancellation frequency at 
nondimensional frequency to = 1.0 (the blade passage frequency) 
for each configuration. Due to the use of different fundamental 
blade passage frequencies as nondimensional reference 
frequencies, mode cut on frequencies are quite different in the two 
cases presented. All mode turn-on frequencies below 
nondimensional frequency to = 5.0 are listed in Table 3.3.2. One 
additional area of difference noted in Table 3.3.1 is the fuselage 
skin thickness; the Gulfstream 695A fuselage skin thickness is 

72 



Table 3.3.1: Long and short wavelength sweep input parameters 
representing the Gulfstream Aerospace Commander 695A and the 
Gulfstream II. Unlisted parameters are as given in Table 2.2.1, 
“Compensated Two Panel Model”. 



Gulfstream 695A 
Long wavelength 
case 

Gulfstream II 
Short wavelength 
case 

propeller RPM 

1518 - 1591 

1305 - 1778 

B 

number of blades 

3 

8 

BPF 

blade passage freauency 

75.9 - 79.6 

174 - 237 
(225 Hz cruise) 

h 

fuselage height 

NA 

7.9' 

d 

(tip to fuselage 
clearance) 

1.18’ 

5.55' 

pnmiin 

0.063" 

0.040" 


8.82’ 

9.00' 

*e 

(external pressure 
propagation 
wavelength) 

19.36' 

8.32' 

Si 

co ref H 

r. 

0.434 

1.257 

S2 

PoH 3 

mi i 

0.085 

0.134 

| S 6 

K e H 

0.3245 

0.7547 

C 

tnanel damning ratios) 

0.01 

0.01 

■■KESMMI 

1.0, 0.6746 

1.0, 0.7263 

Ui 

(nondimensional panel 
velocities) 

solved for 

solved for 

(0 n 

(panel natural 
frequencies) 

0.5, 1.5 

0.5, 1.5 


73 




0.0625", while the Gulfstream II fuselage skin thickness is 0.040". 
The thicker skin of the Gulfstream 695A is 50% more massive than 
the Gulfstream II; in order to keep the simulations representing 
realistic configurations, these data were left as is. However, Figure 
3.3.23 shows the noise reduction predictions for the two aircraft 
with identical skin thicknesses for comparison. 

Table 3.3.2: Mode turn-on frequencies below to = 5.0 for the long 
and short wavelength sweep parameters. 


mode 

# 

Gulfstream 695A 
Long wavelength case 

Gulfstream II 
Short wavelength 
case 

1 

0.905 

0.312 

2 

1.810 

0.625 1 

3 

2.715 

0.937 

4 

3.619 

1.250 

5 

4.524 

1.562 1 

6 


1.874 

7 


' 2.187 j 

8 


2.499 

9 


2.812 

10 


3.124 

1 1 


3.437 

1 2 


3.749 

13 


4.061 

14 


4.374 

1 15 


4.686 

1 6 


i~ 4.999 


Section 3.3.2: Mode Convergence Studxi 

One important decision in numerically solving the governing 
Equations (3.2.22) with its constituent parts is the number of 
acoustic modes required to accurately represent a solution. Use of 


















too few modes will present incorrect results and perhaps hide 
some interesting phenomena. Inclusion of unneeded modes in the 
solution increases computation time unnecessarily. A simple way 
to determine this is to perform a mode convergence and fall-off 
study away from the panel barrier at the ART frequency. Near the 
panel barrier, even those modes which are cutoff will have an 
effect on the resultant pressure (recall that these modes need 
room in the X direction to decay) and as such are important to the 
calculation. Figures 3.3.1 and 3.3.2 show convergence fall off 
results, where the noise reduction is calculated at z = 0.9375 in 
Figure 3.2.1 along the X direction from just inside the panel barrier 
to a point 80 panel lengths downstream, far enough for all cut off 
modes to have decayed. Note that for comparison reasons, all 
falloff results are shown plotted on identical scales. Figure 3.3.1a 
and 3.3.1b show falloff results for the case without external 
pressure sweep. The long wavelength case presented in Figure 
3.3.1a shows a uniformly oscillating noise reduction down the duct. 
For this configuration, only one mode is propagating at the ART 
nondimensional frequency to = 1.0, and this oblique wave bounces 

down the duct, adding and subtracting to the overall noise 
reduction level. Figure 3.3.1b shows the same calculation for the 
short wavelength case, where at the ART frequency three modes 
have turned on and cause a greater fluctuation in the noise 
reduction in the X direction. Figures 3.3.2a and 3.3.2b show the 
same calculation with the external sweep effect, where now the 
noise reduction variation has increased in the long wavelength 
parameters example due to the introduction of more oblique 
waves in the duct caused by the sweeping effect. The calculation 
in the short wavelength case remains qualitatively similar to the 
case without external sweep, except that the variation in noise 
reduction levels is now over a slightly narrower dynamic range. 
Figures 3.3.1 and 3.3.2 indicate that 16 modes generally resolves 
the noise reduction predictions well in the panel near-field for all 




Figure 3.3.1a: Prediction of noise reduction moving 

away from panel barrier for the long wavelength 
case without external sweep. Number of modes in 
calculation as shown in legend. 



Figure 3.3.1b: Prediction of noise reduction moving 

away from panel barrier for the short wavelength 
case without external sweep. Number of modes in 
calculation as shown in legend. 

76 



0 10 20 30 40 50 60 70 80 

Distance from panel barrier (panel lengths) 


Figure 3.3.2a: Prediction of noise reduction moving 

away from panel barrier for the long wavelength 
case with external sweep. Number of modes in 

calculation as shown in l egend. 

2 modes 16 modes 

4 modes 32 modes 



0 10 20 30 40 50 60 70 80 

Distance from panel barrier (panel lengths) 


Figure 3.3.2b: Prediction of noise reduction moving 

away from panel barrier for the short wavelength 
case with external sweep. Number of modes in 
calculation as shown in legend. 
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cases; 32 modes were universally used in all calculations, and 
represents a conservative choice. 

Figure 3.3.3a shows a comparison of the falloff of the 
pressure ratio pt/pE for 32 modes in the duct with external sweep 
in the X direction at height z = 0.9375 for both sets of input data. 

In both cases, the pressure ratio is relatively large near the panel 
barrier due to the effect of cutoff modes. As mentioned 
previously, in the long wavelength case the pressure ratio 
uniformly rises and falls down the duct due to the interaction of 
the one cut on mode with the one dimensional mode. The short 
sweep calculation shows more variation due to the increased 
number of propagating modes. Figure 3.3.3b shows similar 
calculations without the effect of external sweep. Again the long 
sweep case shows a uniform variation as the calculation point 
moves down the duct, and the short wavelength case is less 
regular. Figures 3.3.3a and b predict that the long sweep 
parameters function more effectively to reduce the transmitted 
pressure ratio in the cases with external sweep; without the effect 
of the external sweep, the short wavelength parameters function 
better. 

Section 3.3.3: Comparison of ART anti Identical Tunings 

Figure 3.3.4a shows pressure ratio pt/pE for ART tuned and 
identically tuned panels calculated as a function of frequency at 
location A (x = 2.0, z = 0.9375) in the duct as shown in Figure 3.2.1 
for the long wavelength parameters. This location is two panel 
widths downstream in the X direction, and at a vertical position of 
just under 94% duct height in the Z direction. Identically tuned 
panels have a large transmitted pressure ratio at nondimensional 
frequency = 0.8; this frequency is below the nondimensional 
resonance frequency of 1.0 due to the loading of the apparent 
mass modes and corresponds to a resonance of the combined 
panel-fluid system. At this resonance frequency, the entire panel 
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Figure 3.3.3a: Acoustic pressure ratio p t / pg prediction 

moving away from the panel barrier for the long and 
short wavelength cases with external sweep. 32 
acoustic modes are used in the solution. 



Figure 3.3.3b: Acoustic pressure ratio P t / Pe prediction 

moving away from the panel barrier for the long and 
short wavelength cases without external sweep. 32 
acoustic modes are used in the solution. 
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Figure 3.3.4a: Pressure ratio p t / P E calculation for 

two panel ART tuned and identical panels with long 
wavelength external sweep calculated at Location A 
in Figure 3.2.1. 



Figure 3.3.4b: Pressure ratio p t / p E calculation for 

two panel ART tuned and identical panels with short 
wavelength external sweep calculated at Location A 
in Figure 3.2.1. n n 
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barrier moves in an approximately uniform motion (with the 
exception of the phase variation from panel to panel as a result of 
the external sweep) with large panel velocities. ART panels have 
pressure ratio maximums at nondimensional frequencies of 0.5 
and 1.5 which are much less, since now two resonance frequencies 
are present in the system. Minor perturbations can be seen near 
mode cut on frequencies. Figure 3.3.4b displays results for the 
short wavelength case study; much similarity to the long 
wavelength case is observed with the additional effects of more 

mode cut ons. 

Figures 3.3.5a and 3.3.5b show noise reduction calculations 
for ART and identically tuned panels at the same location in the 
duct for the long and short wavelength case studies, respectively. 
Low frequencies exhibit the usual stiffness dominated behavior, 
and high frequencies show typical mass dominated behavior. As 
can be expected, noise reduction is low at panel natural 
frequencies, and high at the ART design frequency of 1.0. The 
effects of mode cut on can be seen. The long wavelength 
parameters tend to show somewhat better noise reduction 
performance because the ART design frequency is lower and the 
panel sizes are shorter compared to the acoustic forcing 
wavelengths; that is, the longer wavelength case is effectively 

more cutoff. 

Figure 3.3.6a shows the phase difference between panels 3 
and 4 as indicated in Figure 3.2.1 for the long wavelength case.* At 
low frequencies, the phase difference corresponds to that 
associated with the external forcing. (Note that the external 
wavelength of 19.36’ for the long wavelength case divided into 


* Phase difference between panels is calculated from complex velocities 
obtained from the solution of Equations (3.2.22). At each frequency, t e 
sign of the velocity components is checked and a panel phase angle relative 
to the forcing 4> is calculated using the relationship $ = 

arctan(Ui m aginary/Ureal). Successive phase shifts may add multiples of 2* 
to this phase angle difference which do not appear in Figures 3. 3.6a ana 
3.3.6b. 
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Nondimcnsional Frequency 

Figure 3.3.5a: Noise reduction calculation for two panel 

ART tuned and identical panels with long wavelength 
external sweep calculated at Location A in Figure 3.2.1. 



Nondimensional Frequency 

Figure 3.3.5b: Noise reduction calculation for two panel 

ART tuned and identical panels with short wavelength 
external sweep calculated at Location A in Figure 3.2.1. 
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Nondimensional Frequency 

Figure 3.3.6a: Phase difference between panels 3 and 4 

as a function of nondimensional frequency for the 
long sweep wavelength case. 



Nondimensional Frequency 

Figure 3.3.6b: Phase difference between panels 3 and 4 

as a function of nondimensional frequency for the 
short sweep wavelength case. 
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360 degrees yields about 19 degrees phase shift between panels; 
this is the low and high frequency phase difference limit.) Mode 
cut on causes one irregularity in the phase difference prediction at 
co ~ 0.9; the other anomaly occurs near the high panel resonance 
frequency. Figure 3.3.6b shows a similar calculation for the short 
wavelength case; here additional mode cut ons (3 modes cut on 
between the panel resonance frequencies of co = 0.5 and 1.5) 
contaminate the phase difference calculation. Low and high 
frequency phase difference limits have changed to a larger limiting 
value due to the shorter external forcing wavelength. Figure 
3.3.7a and 3.3.7b show panel velocity relationships for the long 
and short wavelength ART cases. Both plots indicate that the panel 
systems have matching velocity amplitudes at the ART design 
frequencies co = 1.0. The long wavelength case indicates that the 
panels reach higher velocity amplitudes around the natural 
frequencies than the short wavelength case. 

Figures 3.3.8a and 3.3.8b show a comparison of ART tuning 
and the average noninteracting noise reduction results (ANNR). 
Recall that the ANNR is a idealized measure of the average noise 
reduction present in the system where only uniform panels are 
used, and the resulting acoustic fields do not interact, as shown in 
Chapter 2, Figure 2.4.5. This comparison allows the mass law and 
stiffness law effects associated with moving the panel resonance 
frequencies apart to be separated from the effects of acoustic 
cancellation or cutoff which are critical to optimum ART 
performance. In both the long and short wavelength cases, ART 
tuning performs better than the ANNR calculation. In the long 
wavelength case, between 10 dB and 25 dB of increased noise 
reduction is achieved. The short wavelength case, while still 
outperforming the ANNR result, is not as dramatic as the long 
wavelength case, providing an 8 dB to 12 dB noise reduction 
increase. 

As a final comparison of ART and identical tunings. Figures 
3.3.9a and 3.3.9b show a comparison of the individual pressure 
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Figure 3.3.7a: Velocity relationships for panels 3 and 4 

as a function of nondimensional frequency for the 
long sweep wavelength case. 



Figure 3.3.7b: Velocity relationships for panels 3 and 4 

as a function of nondimensional frequency for the 
short sweep wavelength case. 
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Nondimensional Frequency 

Figure 3.3.8a: Averaged noninteracting noise reduction 

for identical panels and two panel ART noise reduction 
prediction for long sweep wavelength case calculated at 
Location A in Figure 3.2.1. 



Figure 3.3.8b: Averaged noninteracting noise reduction 

for identical panels and two panel ART noise reduction 
prediction for short sweep wavelength case calculated at 
Location A in Figure 3.2.1. 
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Figure 3.3.9a: Relative comparison of the first 21 

modal contributions to noise reduction in the 
long sweep wavelength case at the ART design 
frequency. Note that the comparison is logarithmic. 
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Figure 3.3.9b: Relative comparison of the first 21 

modal contributions to noise reduction in the 
short sweep wavelength case at the ART design 
frequency. Note that the comparison is logarithmic. 
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mode coefficients for the calculation of the acoustic pressure in the 
duct at the ART nondimensional design frequency. Since there is a 
wide dynamic range of these pressure coefficients and because 
there is a net noise reduction across the panel barrier, the negative 
logarithm of the pressure coefficient results is taken. This 
mathematical interpretation means larger values of the modal 
pressure coefficients contribute more to noise reduction in the 
system; that is, less transmitted noise through the panel barrier 
occurs with larger values of the modal pressure coefficients. The 
following comments are made in that spirit. While the ART tuning 
modal pressure coefficients are not always greater than the 
uniform panel tuning coefficients over the range of mode numbers 
shown, the ART results for lower propagating modes are always 
greater than similar results for identical panels. Note that for each 
propagating pressure mode (as shown in Figure 3.3.9a, modes 0 
and 1 for the long sweep wavelength case, and modes 0 through 3 
for the short wavelength case shown in Figure 3.3.9b), the ART 
tuning methodology contributes to noise reduction; in some cases, 
this contribution is more than one order of magnitude. Also note 
that modes 8 and 16 have large values because these modes are 
integer multiples of the number of panels in the system. In 
general, as the mode number increases, the negative logarithm of 
the pressure coefficient increases more for ART panels than for 
identical panels, indicating that ART redirects energy into the 
higher, cutoff modes. 


88 



Section 3.3.4; Model Parameter Studi&a 


A comprehensive set of parameter studies provides further 
insight with respect to tuning an ART panel array in the presence 
of an external propagating pressure field. Note that for Figures 
3.3.10 through 3.3.14, noise reductions are calculated at position A 
as shown in Figure 3.2.1. 

Figure 3.3.10a shows noise reduction calculations in the duct 
for three separate sets of input parameters, where the higher 
panel natural frequency is increased in order to increase the noise 
reduction bandwidth. In the long wavelength case, this is quite 
beneficial, since mode turn -on frequencies do not occur as often as 
the short wavelength case. Splitting the natural frequencies not 
only increases the noise reduction bandwidth, but it also increases 
the peak noise reduction obtained. Figure 3.3.10b shows that this 
procedure is still beneficial with the short sweep wavelength 
parameters, even in the presence of many mode turn-ons. This 
study suggests that when tuning an ART panel array, it is most 
desirable to attempt to keep mode turn-on frequencies away from 
the desired frequency of maximum noise reduction; however, this 
may not always be feasible. In general, increasing the ART noise 
reduction bandwidth by separating panel natural frequencies up to 
a difference of nondimensional frequency of 2.0 is most always 
beneficial. (Chapter 5 will show experimental verification of this 
fact.) 

Figures 3.3.11a and 3.3.11b show the effect of increased 
panel damping on the model. As with the four panel model 
presented in Chapter 2, increased damping tends to increase noise 
reduction levels at the panel natural frequencies, as each panel in 
effect becomes less acoustically transparent. Since the panels are 
operating in a regime father away from perfect out-of-phase 
behavior, absolute maximum noise reductions are less. However, 
the general frequency bandwidth at a given noise reduction 
amount remains the same, indicating the robustness of the ART 

89 




Figure 3.3.10a: Parameter studies with two panel tuning 

varying the higher natural frequency as indicated above 
to increase the noise reduction bandwidth in the long 
wavelength case. 



Figure 3.3.10b: Parameter studies with two panel tuning 

varying the higher natural frequency as indicated above 
to increase the noise reduction bandwidth in the short 
wavelength case. 
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Nondimensional Frequency 


Figure 3.3.11a: Parameter studies with the two panel model 

varying the panel damping ratio as indicated above in the 
long wavelength sweep case. 



Nondimensional Frequency 


Figure 3.3.11b: Parameter studies with the two panel model 

varying the panel damping ratio as indicated above in the 
short wavelength sweep case. 
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tuning concept, and that ART tuning might be incorporated into 
other modes of noise reduction such as increased damping. 
Furthermore, in the short sweep wavelength parameters case, 
increased damping seems to play an even lesser role, indicating 
the tendency of the external forcing to dominate system behavior. 

Figure 3.3.12a and 3.3.12b show the effect of the vanation of 
the apparent mass loading on the ART panel array. In both the 
long and short wavelength cases, increasing the apparent mass 
loading tends to couple the panels together through the 
surrounding fluid such that approximate out-of-phase behavior is 
no longer possible. As such, maximum noise reduction and noise 
reduction bandwidth are simultaneously reduced. This effect is 
particularly noticeable in the long sweep wavelength case study, 
where over 25 dB of noise reduction is lost between normal air 
mass loading and eight times normal air mass loading. If operation 
of an ART panel array were necessary in a heavy fluid loading, it 
would become necessary to design for this loading. An alternative 
option might be to increase panel mass in order to maintain the 

apparent mass to panel mass ratio. 

Figures 3.3.13a and 3.3.13b show the noise reduction 
predictions obtained by variation of the parameter A as denoted 
by Equation (3.2.17) in Section 3.2. Recall that when A — 1, fluid 
loading on only one panel side is considered, whereas when A - 2, 
fluid loading on both sides on the panel is considered. As 
expected, the solutions are quite similar except near panel natural 
frequencies. Panel velocity is the greatest around the panel 
natural frequency. Variation of A as shown effectively doubles the 
apparent mass loading, and when the panel velocities are large, 
this loading is most critical. Therefore, for all calculations, 
apparent mass loadings on both panel sides are routinely 
considered. 

Propeller rotation direction can be controlled by the sign of 
the nondimensional parameter KeH (S6 in Table 3.3.1), the default 
value used herein has been an upward or positive propagation 
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Nondimensional Frequency 


Figure 3.3.12a: Parameter studies with the two panel model 

varying the panel apparent mass loading as indicated above 
in the long wavelength case. 



, Nondimensional Frequency 

Figure 3.3.12b: Parameter studies with the two panel model 

varying the panel apparent mass loading as indicated above 
in the short wavelength case. 
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Figure 3.3.13a: Noise reduction calculation for two panel 

ART tuned array and long wavelength sweep parameters 
with A = 1 and A = 2. Noticeable difference occurs only at 
resonance frequencies. 



Figure 3.3.13b: Noise reduction calculation for two panel 

ART tuned array and short wavelength sweep parameters 
with A = 1 and A = 2. Noticeable difference occurs only at 
resonance frequencies. 
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direction with respect to the coordinate system shown in Figure 
3.2.1. Reversal of the external sweep can be accomplished by use 
of negative values for the parameter K e H. Figures 3.3.14a and 
3.3.14b shown the dramatic effect that the propagation direction of 
the external forcing can have upon the noise reduction across the 
panel barrier with respect to a given location; again, noise 
reductions are calculated at Location A as shown in Figure 3.2.1. 

For the latter location, reverse sweep shows significant increases m 
noise reductions at certain frequencies. Using the long wavelength 
sweep parameters, as shown in Figure 3.3.14a, noise reduction is 
increased at higher frequencies between 3 and 8 dB. Also note 
slightly better noise reductions between the panel resonance 
frequencies. Similar results are seen in Figure 3.3.14b, 
representing the short external wavelength parameters. The 
difference in these cases might be attributed to different 
interactions of the panel barrier with the upper and lower duct 
boundaries. In the upper sweep propagation mode, the exterior 
propagating pressure last sees panel number 8 (referring to Figure 
3.2.1), which is ART tuned to a higher natural frequency and has 
less mass. In the downward sweep direction, the propagating 
pressure last sees panel number 1, ART tuned to a low natural 
frequency with more mass. The different sweep directions set up 
different groups of oblique waves which propagate down the duct, 
and the increase in noise reduction can be attributed to a greater 
local cancellation effect. Noise reduction calculations at a point 
symmetrically opposite of Location A (x = 2 . 0 , z = 0 . 0625 ) appear 
worse. This will be verified and discussed in the next section. 

Section 3.3.5: Pressure and Noise Rfflllftiftl) Mappi ng 

The solution of the governing Equations (3.2.22) yields a full 
acoustic solution; pressures (and therefore noise reduction levels) 
can be calculated at all locations within the duct knowing the panel 
velocities. Pressure and noise reduction mappings at one value of 
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Nondimensional Frequency 


Figure 3.3.14a: Forward and reverse sweep direction noise 

reduction calculation for two panel ART tuning with long 
wavelength external sweep calculated at Location A. 



Nondimensional Frequency 


Figure 3.3.14b: Forward and reverse sweep direction noise 

reduction calculation for two panel ART tuning with short 
wavelength external sweep calculated at Location A. 
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frequency throughout the entire duct domain indicate the sound 
levels to be expected at all locations. Such calculations indicate the 
nature of the sound fields within the duct at a desired frequency, 
specifically, “hotspots” and quiet areas can be identified.. 

Figure 3.3.15a shows the pressure mapping in the duct at the 
ART frequency co = 1.0 for the long sweep wavelength parameters.. 
Note that the three dimensional transmitted pressure ratio plot has 
X location in units as nondimensionalized by panel height H, Z 
location in the duct (vertical position per Figure 3.2.1) on the “y” 
axis (0 < z < 1), and acoustic pressure on the “z” axis. Rather high 
pressure ratio values are predicted in the upper portion of the 
duct closest to the panel barrier in the long sweep wavelength 
case. The pressure ratio tends to decay to a somewhat uniform 
value after movement away from the panel barrier of 2-3 panel 
widths. The short wavelength parameters pressure mapping in 
Figure 3.3.15b shows the effect of more propagating modes at the 
ART frequency co = 1.0, where three distinct pressure peaks are 
observed due to the three propagating modes. Again the 
uppermost portion of the duct has higher transmitted pressure 
ratios, where noisy regimes are to be expected. 

Sound mapping in decibels provides a more intuitive 
interpretation of the relative sound levels in the duct. Figure 
3.3.16a plots noise reduction in the duct for the long sweep 
wavelength parameters case. When translated into decibels, the 
ART effect is quite dramatic; after moving a distance of 5 panel 
diameters away from the barrier, at least 20 dB of noise reduction 
has occurred. Also, the presence of only one propagating mode 
insures a more uniform noise reduction throughout all areas of the 
duct. Figure 3.3.16b predicts that the short wavelength 
parameters produce more variation in the noise reduction 
throughout the duct. Figures 3.3.16c and 3.3. 16d show the same 
calculations without the effect of external sweep. In the long 
sweep wavelength case in Figure 3.3.16c, the noise reduction has 
approached a “tent-like” structure, appearing relatively regular in 
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Figure 3.3.16a: Noise reduction map for the ART tuned panel 

array at the maximum ART noise reduction frequency co - 
in the long sweep wavelength case. 




Figure 3.3.16b: Noise reduction map for the ART tuned panel 

array at the maximum ART noise reduction frequency co = 1.03 
in the short sweep wavelength case. 
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Figure 3.3.16c: Noise reduction map for the ART tuned panel 

array at the maximum ART noise reduction frequency to = 1.0 
for the long wavelength parameters without external sweep. 
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Figure 3.3. 16d: Noise reduction map for the ART tuned panel 

array at the maximum ART noise reduction frequency fl> — 1.03 
for the short wavelength parameters without external sweep. 
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construction. Figure 3.3. 16d suggests again that the short 
wavelength parameters produce a wider variation in noise 
reduction throughout the duct. Some of the propagating pressures 
can interact to produce quite dramatic, but only locally high 
effectiveness in attenuating the transmitted noise. Nonetheless, 
the ART treatment does produce a substantially quieter interior 

noise level. 

Figures 3.3.17a through 3.3.17d show cross sectional noise 
reduction calculations at the indicated panel widths downstream 
from the panel barrier calculated at to = ©ART- Note in these figures 
that noise reduction is now placed on the abscissa, and vertical 
location in the duct is plotted on the ordinate. In Figure 3.3.17a, 
the long sweep wavelength parameters produce a more regular 
noise reduction level moving further away from the wall. This is 
in contrast to the short sweep wavelength parameters shown in 
Figure 3.3.17b, which produce large variations across the duct. 
Figures 3.3.17c and 3.3. 17d show the equivalent cases without the 
sweeping effect; in the long sweep wavelength parameters, the 
“tent-like” noise reduction prediction of Figure 3.3.16c is examined 
Note the decay of pressure (noise reduction increase) indicated on 
the centerline as the calculation point moves away from the panel 
barrier. Also, the noise reduction calculation is not symmetrical; 
there is variation in the upper and lower portions of the duct due 
to the different boundary condition represented by panels of 
varying mass and natural frequency. 
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Z Location 



Noise Reduction (dB) 

Figure 3.3.17a: Cross sectional noise reduction predictions 

at indicated X locations for the long wavelength sweep case 



Noise Reduction (dB) 

Figure 3.3.17b: Cross sectional noise reduction predictions 

at indicated X locations for the short wavelength sweep case 
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Z Location ^ Location 



Noise Reduction (dB) 

Figure 3.3.17c: Cross sectional noise reduction predictions 

at indicated X locations for the long wavelength sweep 
parameters without external sweep. 



Noise Reduction (dB) 

Figure 3.3. 17d: Cross sectional noise reduction predictions 

at indicated X locations for the short wavelength sweep 
parameters without external sweep. 
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Section 3.3.6; Comparison to the Two Panel Mo tel 


Having derived a complete acoustic solution for the model 
presented in this chapter, a comparison to the simplified model of 
Chapter 2 provides additional insight into the model limitations 
and strengths. Figures 3.3.18a and 3.3.18b show the Chapter 3 
model reduced to a panel barrier containing only two panels with 
no exterior pressure sweeping effect, and comparisons are made in 
the far field (5 duct widths away) and near field (1 panel width, or 
a half duct width away). The far field parameters used in the 
Chapter 3 model in Figure 3.3.18a shows excellent agreement with 
the Chapter 2 model, except around nondimensional frequency co = 
3.0 where the first mode turn on occurs for the 2 panel geometry. 
Figure 3.3.18b calculated in the near field shows similar 
disagreement at the mode turn-on frequency, except that the 
variation is larger since the calculation point is much closer to the 
panel barrier. These variances between the Chapter 2 model 
(which assumed all higher modes decayed and only the acoustic 
zero mode propagated) and the Chapter 3 model are to be 
expected. Note that the maximum noise reduction prediction for 
both models in the near and far fields is quite comparable. 

Another interesting case is the comparison of the two panel 
Chapter 2 model with the full 8 panel Chapter 3 model in the near 
and far fields, as shown in Figure 3.3.19. Including a panel barrier 
with more panels creates more mode turn ons, which lower the 
noise reduction predictions at the ART design frequency and 
causes more subtle variations in the noise reduction prediction as 
frequency increases due to higher mode turn ons. 

Section 3.3.7; Miscellaneous — Resul ts 

A number of additional effects may be explored with the full 
acoustic solution for the external pressure sweeping model. Figures 
3.3.20a and 3.3.20b show noise reduction results for a high value 
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Figure 3.3.18a: Far field (5 duct widths) noise reduction 

predictions for the Chapter 2 model (with decay of 
higher modes) plotted with the Chapter 3 model including 
all modal contributions for the 2 panel configuration. 



Figure 3.3.18b: Near field (1 panel width) noise reduction 

predictions for the Chapter 2 model (with decay of 
higher modes) plotted with the Chapter 3 model including 
all modal contributions for the 2 panel configuration. 
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Nondimensional Frequency 

Figure 3.3.19a: Far field (5 duct widths) noise reduction 

predictions for the Chapter 2 model (with decay of 
higher modes) plotted with the Chapter 3 model including 
all modal contributions for the 8 panel configuration. 



Nondimensional Frequency 

Figure 3.3.19b: Near field (1 panel width) noise reduction 

predictions for the Chapter 2 model (with decay of 
higher modes) plotted with the Chapter 3 model including 
all modal contributions for the 8 panel configuration. 
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Nondimensional Frequency 


Figure 3.3.20a: High damping ratio noise reduction 

prediction for the long wavelength sweep case 
compared to the averaged noninteracting noise 
reduction prediction. 
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Nondimensional Frequency 


Figure 3.3.20b: High damping ratio noise reduction 

prediction for the short wavelength sweep case 
compared to the averaged noninteracting noise 
reduction prediction. 
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of panel damping (£ = 0.4) compared to the ANNR results for 
identical panels. Long wavelength parameter results are shown in 
Figure 3.3.20a; short wavelength parameter results are shown in 
Figure 3.3.20b. In both cases, ART performs about 5 to 7 dB better 
than uniform panels in the design frequency range around 
nondimensional frequency = 1.0. This again suggests that the ART 
concept might be implemented with more traditional noise 
reduction techniques such as increasing panel damping to achieve 
better quieting performance. 

In practice, it is impossible to have exactly matching panel 
properties. Panel natural frequencies and damping ratios will 
always show some performance variations, although perhaps over 
a small margin. Figures 3.3.21a and 3.3.21b show the effect of this 
“detuning” on ART noise reduction, where panel mass ratios, 
natural frequencies, and damping ratios have been varied 
randomly by ±10% as shown in Table 3.3.3. In both the long and 
short wavelength cases, the ART tuning methodology is quite 
robust, and these tuning variations cause no negative effects. 

There are some additional ART noise cancellation peaks seen in 
Figures 3.3.21a and 3.3.21b due to the ART effect occurring 
between detuned parameters. For example, in Figure 3.3.21a, 
there is an ART peak occurring at a nondimensional frequency just 
above the low natural frequency panel resonance minimum. This 
ART peak is a result of the interaction of the panels with resonance 
frequencies of 0.45 and 0.55 respectively. In all other respects, 
the ART noise reduction performance is essentially the same as the 
cases where only uniform panels are used. 




Figure 3.3.21a: Noise reduction calculation at Location A 

for detuned panel properties in the long wavelength sweep 

case. 



Nondimensional Frequency 

Figure 3.3.21b: Noise reduction calculation at Location A 

for detuned panel properties in the short wavelength sweep 

case. 
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Table 3.3.3: Panel detuned parameter values used in Figure 3.3.22. 


panel # 

MRj 

COn 

Ci 

1 

0.9 

0.55 

0.01 1 

2 

0.603 

1.35 1 

0.011 

3 

1.1 

0.45 

0.009 

4 

0.737 

1.65 

0.009 

5 

1.1 

0.55 

0.011 

6 

0.737 

1.35 

0.009 

7 

0.9 

0.45 

0.011 

8 

0.603 

1.65 

0.009 


“Off-tuning” considerations are another topic which may be 
explored with the Chapter 3 model. Recall that any number of 
panels may be used in the model, and these panels may be tuned 
in any manner. Occasions may arise in real structures where it 
will be impossible to have an even number of panels, and hence 
complete ART panel pairs. Figure 3.3.22 shows noise reduction 
predictions for a barrier with 9 panels comprising four ART tuned 
pairs, and one extra panel with a natural frequency of co = 0.5 
placed at the uppermost portion of the panel barrier. Sweep 
propagation remains in the +z direction. Figure 3.3.22a shows the 
noise reduction predictions for the long sweep wavelength 
parameters, and two cases are shown at the locations indicated in 
the legend. These positions correspond to locations which are 
symmetry opposite the duct centerline at z = 0.5. Note that the 
noise reduction prediction at the lower position (x = 2.0, z — 0.06) 
shows an additional noise reduction peak. Higher frequency noise 
reduction is also enhanced at the lower location. Similar results 
may be seen in Figure 3.3.22b, where the short wavelength 
parameters have been employed. Figures 3.3.22c and 3.3.22d 
show cross sectional noise reduction predictions at coart — 1*0 for 
the long and short wavelength parameters cases. These figures 
may be directly compared to Figures 3.3.17c and 3.3.16d, which 
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Nondimensional Frequency 

Figure 3.3.22a: Modelling a panel wall irregularity for the 

long sweep wavelength case. 
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Nondimensional Frequency 

Figure 3.3.22b: Modelling a panel wall irregularity for the 

short sweep wavelength case. 
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Figure 3.3.22c: Cross sectional noise reduction predictions 

at indicated X locations for the long wavelength sweep case 
with panel tuning irregularities. Calculation atw^ T -1.0. 



Noise Reduction (dB) 

Figure 3.3. 22d: Cross sectional noise reduction predictions 

at indicated X locations for the short wavelength sweep case 
with panel tuning irregularities. Calculation at(o ART = 1.0. 
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show the same calculation for an even number of ART tuned pairs. 

In the long wavelength case (Figures 3.3.17c and Figure 3.3.22c), 
the use of an extra panel has increased the minimum noise 
reduction found in the duct by about 5 dB, and increased the 
maximum noise reduction by 2 dB. Slightly more impressive 
results may be seen in the short wavelength parameters case 
(Figures 3.3. 17d and 3.3.22d). In the latter example, the noise 
reduction calculation at 4 panel widths downstream has increased 
by about 10 dB at some locations, although absolute maximum 
noise reductions have decreased slightly. This demonstrates the 
complex interaction of the panels in the panel barrier and the 
sound fields which are produced as a result of the panel motions, 
and implies that some noise reduction gain may be achievable by 
considering “off-tuning” possibilities. 

Recall that the long wavelength parameters (from the 
aircraft quoted by Wilby and Wilby 24 ) are obtained based on a 
skin thickness of 0.0625". The short wavelength parameters 
(obtained from a study by Kuntz and Prydz 25 ) assume a skin 
thickness of 0.040". Figure 3.3.23 is used to compare the noise 
reduction predictions for the long and short sweep wavelength 
parameters with the two different apparent mass to panel mass 
ratios. Mass domination is a high frequency concept due to the 
imco component of the panel impedance formula expressed by 
Equation (2.2.3). As shown in Figure 3.3.23a and 3.3.23b, this 
extra mass is responsible for about a 5 dB increase in noise 
reduction at high frequencies, and about a 2-3 dB difference at the 
ART design frequency. 

Four panel tuning is also a possibility for the model 
developed in this chapter. Figure 3.3.24a shows noise reduction 
predictions for four ART-tuned panels compared to identical 
panels. This long wavelength four panel ART tuning sweep case 
predicts increases in noise reduction of about 35, 5, and 5 dB at the 
three ART design frequencies of to = 1.0, 2.0, and 3.0. This four 
panel tuning is not as successful in the case with the short 
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Nondimensional Frequency 

Figure 3.3.23a: Comparison of noise reduction predictions 

with two apparent mass loading values for the long 
wavelength parameters sweep case. 



Nondimensional Frequency 


Figure 3.3.23b: Comparison of noise reduction predictions 

with two apparent mass loading values for the short 
wavelength parameters sweep case. 
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Nondimcnsional Frequency 

Figure 3.3.24a: Noise reduction calculation for four panel 

ART tuned and identical panels with long wavelength 
external sweep calculated at Location A in Figure 3.2.1. 



Nondimensional Frequency 

Figure 3.3.24b: Noise reduction calculation for four panel 

ART tuned and identical panels with short wavelength 
external sweep calculated at Location A in Figure 3.2.1. 
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wavelength parameters, however, as shown in Figure 3.3.24b, 
where the increases are about 25, 4, and 1 dB respectively. Recall 
that the short wavelength parameters are operating in a higher 
frequency regime, and as such, are less effective. Figures 3.3.24a 
and 3.3.24b show that while the ART concept can produce N-l 
noise cancellation effects for N different panel tunings in the near 
field, it is inherently more effective at lower frequencies where 
the panel systems are more acoustically compact; as such, the two 
panel tuning is more worthwhile. Additionally, existing noise 
reduction methods tend to work better at higher frequencies and 
are less effective at lower frequencies. It is in the latter frequency 
domain where the employment of the ART technique can achieve 
the most dramatic results. 


119 



Chapter 4 

The General ABT Pangl Analysis 


Section 4.1: Introduction 

In the analyses presented thus far, panel impedances have 
been represented by the familiar relation 

Z ra = Ri + i|Mi(0 + 5i) (4.1.1) 

This mathematical statement of the panel dynamics assumes that 
the panel itself is represented by a rigid flat plate section; no 
higher structural modes are represented here. While this 
representation is analytically much easier to deal with, real panels 
attached to a stringer-rib frame are generally held in place using 
rivets, and this riveting creates a panel mode shape which is 
“bowed” to some degree depending upon the actual boundary 
condition. It is necessary to investigate the effect such panel mode 
shapes might have on ART performance. One particularly 
interesting possibility is the contribution of higher structural 
modes to the ART process. 

Section 4.2: Analysis — Derivation 

The most general analytic configuration is shown in Figure 
4.2.1. This geometry simply includes a standard four panel array. 
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1 1 

12 ^ 

21 

22 


Any desired 
structural shape, 
resulting from 
any desired 
boundary 
condition. 


Figure 4.2.1: Configuration for ART 

Real Panel Analysis 


two panels on top and two panels on the bottom. Four panels for 
each identical panel subsystem have been chosen as a convenient 
number for equation derivation, but in reality, any rectangular 
array of panels may be considered. For the present case, panels 
are numbered in the standard matrix-like notation as shown. Note 
that no restriction has been placed on the boundary conditions of 
each panel. The latter assumption allows the analysis to proceed 
most generally, and at a later appropriate time, the boundary 
conditions will be chosen through specification of the appropriate 
structural mode shapes for the panels. The panel array is 
considered to be one subsystem mounted on an infinite wall 
composed of infinitely many panel subsystems. An anechoic 
termination is initially assumed beyond the wall. 
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To formulate the governing equations, consider panel 11 
(note that this is an arbitrary choice). A general acoustic velocity 
function may be defined as 


u(x,y,z,t) = XlUjj \|/ij(y,z) e i(<ot ' kxi i x> (4.2.1) 

i j 

where the subscript indices i and j refer to acoustic modes, and 
yij(y,z) describes the spatial dependence of an acoustic modal 
function. Ujj are the corresponding complex acoustic modal 
velocity amplitudes. Note that the acoustic modes are assumed 
harmonic in time and propagating in the +x direction. Similarly, a 
general acoustic pressure may be expressed as 

p(x,y,z,t) = X X Pij Vij(y.z) e i(<ot ' kxi * x) (4.2.2) 

i j 

with Pij representing the complex pressure amplitude of acoustic 
mode ij. Note that the relationship 

Pij = Po _a_ Ujj (4.2.3) 

kxij 

may be defined as a consequence of the x-component of the 
momentum equation. 

For the structural problem, the wall displacement may be 
denoted as 


Wy*t) = II A{,‘ OpVy-z) e M (4.2.4) 

p q 

where Apq 11 is the complex structural modal amplitude for mode 
pq of panel 1 1 and <|>pq ^ is the individual panel shape function for 
mode pq of panel 11. A convention is made here that superscripts 
will always refer to the plate location as per Figure 4.2.1; 
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subscripts will always refer to either acoustic modes using i and j, 
or structural modes using p and q. Note that the overall wall 
shape function y and the individual panel shape function <t> are 

separable functions; that is, for panel 11, 

<t>pVy> z ) = 

vW y,z) = ^p(y) ^q( z ) 


Differentiation of Equation (4.2.4) results in an expression for the 
wall velocity as 

u w = ^y,z,t) = 1 1 ioApq 4'^y.z) «>“ (4-2.5) 

at p q 

The acoustic and structural equations are linked at the 
panel/fluid interface, or at x = 0, where the acoustic velocity is 
equivalent to the wall velocity 

“acoustic^-yAO = U wall (0,y,z,0 (4.2.6) 

This boundary condition encourages placement of Equation (4.2.1) 
equal to Equation (4.2.5) at x = 0, resulting in 

IlUjj Vij(y,z) = XXio>ApJ<|>pq(y,z) (4.2.7) 

i j V q 

The acoustic modes and the structural modes each form an 
orthogonal set with respect to themselves. However, the acoustic 
modes are not necessarily orthogonal to the structural modes. 
Application of orthogonality to Equation (4.2.7), and substituting 
dS = dydx yields 
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X X Uij 
i j 



Vij Vij dS total = X S iwA 



Vrs <t>p 1 c flSii( 4 . 2 . 8 ) 


The only modes of the left hand side of Equation (4.2.8) which 
survive are ^rs- Therefore, with a change in indices, 

XXicoAjjjf Vij^pqdSn 

TT n p q hn ( 4 . 2 . 9 ) 

u ‘j f 

I \|/jj 2 dStotal 

J$ total 


By virtue of the momentum equation invoked to produce Equation 
(4.2.3), the complex acoustic pressure amplitudes may be 
expressed as 


IIipio 2 f Vij Pp' q d Sn 


P <1 Jsil AH 

(4.2.10) 

— — - i-vpq 

kxij I Vij 2 dStotal 

/stotil 



A convenient notational shorthand allows the expression of 
Equation (4.2.10) in a more compact form 


where 



= X STjjpq Apq 

p q 


(4.2.11) 



1 

•o’ 

e 

to 

[ Vij ^pq dS l 1 


Js\\ 

^Xij 

Xj/jj 2 dStotal 

4 

StoUl 


(4.2.12) 
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Note an additional notational convention where S p u is the area of 
the single panel 11; its associated derivative dS p u in integral 
equations implies integration over the surface of panel 11. 
Similarly, S to tal is the area of the entire panel array (equal to 4 
times the panel width W times the panel height H); the associated 
derivative dS to tal in integral equations implies integration over the 
entire panel assembly. 

Lagrange's equations may be used to derive the equations of 
motion for the panel structure itself. 26 In general 

Mpq [(A P q e 1C0t )u + 2CpqC0^ p 1 q {A P q e ,wt )t] (4.2.13) 

+ =QiV i “ 

where Mpq 11 is the generalized mass of mode pq of panel 11, Cpq 11 
is the damping ratio of mode pq of panel 11, Q>npq^ is the 
undamped natural frequency of mode pq of panel 11, and Q P q !1 is 
the sum total of all generalized forces acting on panel 11. 

In particular, the generalized forces are worth further 
perusal. The right hand side of Equation (4.2.13) without the 
harmonic excitation e io)t may be expressed in general for any 
panel ap as 


«/ P p q = ss<“ b +^ < 4 - 2 - 14 > 

a b 

where the term Qe is the external forcing term, resulting from 
holding the wall fixed and observing the incident wave exhibit a 
hard wall reflection on panel ap. The summation term is a 
convenient way of representing the generalized force contribution 
of all panels in the entire panelled system on panel aP where the 

analysis is focused. The notation represents the effect of the 
motion of mode pq of panel ab on panel ap. In more general form, 

Equation (4.2.14) may be written for the four panel analysis 
configuration focusing on panel aP as 
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(4.2.15) 



Here the subscript m indicates a generalized force resulting from 
motion of panel ab on <xp. Equations (4.2.14) and (4.2.15) are a 
direct result of the concept of superposition applied to the panel 
subject to analysis. 

Note that the generalized force on the right hand side of 
Equation (4.2.15) may be written as 



pab(0,y,z) may be replaced using Equation (4.2.2) as 


p ab (0,y,z) = 11??/ ¥,/y.z)dSp b < 4 - 2 - 17 > 

» j 

Substitution of Equation (4.2.11) for Pjj a b and back substitution of 
Equation (4.2.17) into (4.2.16) yields the generalized force on panel 
<xp as a result of the motion of panel ab as 


= ( <$(y,z) XLIlO A*5 Vij(y.z) isf (4.2.18) 

] i j c d 

Furthermore, substitution into Equation (4.2.18) yields 

[ »p 5 1 1 ¥ij dS“ P = [ V?j dS Twa i * iTfp, (4- 2 - 1 9 > 

J i j IPO <» 2 1 

where a panel independent quantity Jij may be written as 
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(4.2.20) 


Jij = 2ESL ( <(ljj dStotal 
ip 0 co 2 ] 

Therefore, the entire generalized force on aP will consist of one 
contribution from each panel, and may be most compactly 
represented as 


sscC b+ <=E2: 

a b a b 


Lc d i j 


IJL'U 


(4.2.21) 


With the help of Equation (4.2.13) expressed in the most 
general form for any panel, the panel governing equations may be 
written as a linear system of the form 


XXXX 

a p p q 


^pq IYI pq 


-XX 

a b 




c d 


-o“ p 


(4.2.22) 


where M a P pq is the result of using Lagrange's equations to 
represent the panel dynamics. 



H- 


+ CD 


f 

Ilpq 


(4.2.23) 


QE a ^pq is the external forcing, and the complex modal panel 
amplitudes A a Ppq are the solutions to the linear system. 

Expression of Equation (4.2.22) in nondimensional terms 
allows for the most general engineering interpretations. Note that 
in the combination of Jij, T“Pijpq, and T ab ij C d. some simplification 
occurs; also, the external generalized force is now stated explicitly 

in nondimensional terms, resulting in 
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yp zp sm sm 

SEES 

a=l p.l p=l q=l 


yp zp sm sm 

'ctppq ' X X X X Tlabcdappq = K appq| 

a=l b= 1 c=l d=l 

(4.2.24) 


where the nondimensional terms v,T], and k are denoted by 


v -l Sl lA ap M ap 
v aPpq ~ A pq ^Pq 


< +2iC?5 CO CO 2 


(4.2.25) 


.-27 ab 3 ” 1 a ™7«P , , Q «P / 0) r \ 
Tlabcdappq “ 1C0 A cd ^PQ Vij dS p ^ j 


I <t>p S q Vij dS a p b 




ii dStotal 

(4.2.26) 


K aPpq — 


i 


”~aB 

<b dS 
Y pq 


ap 

P 


(4.2.27) 


Additionally, from Equation (4.2.24), yp is the number of panels in 
the model in the y direction (assumed to be equal to 2), zp is the 
number of panels in the z direction (assumed to be 2), and sm is 
the number of structural modes to be included in the calculation. 
From Equation (4.2.25), other nondimensional variables include 


Si = 


c 


C _pWH 2 

02 = — — 

mu 

where a> r is a user-chosen reference frequency, usually set to the 
lowest ART cancellation frequency, H is the panel height, c is the 
speed of sound in air, p is the density of air surrounding the 
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panels, W is panel width, and mi 1 is the mass of panel 1 1 (a 
reference panel mass) as shown in Figure 4.2.1; 



0>r 


A «P 

A pq P P C 


Pe 


is the complex nondimensional modal amplitude for mode pq of 
panel ap, where Pe is the total external pressure on panel ab equal 
to (2Pi - Pt). Pi is the incident pressure on the upstream panel 
side, and P t is the transmitted pressure through the panel array. 
Other nondimensional relationships from Equation (4.2.25) are 
defined as follows: 


Mpq = 

1 1 
1 1 


m 


ap 


m 


< 


CO 


ap 

n pq 


COr 


©■-«- 

C0 r 


From Equation (4.2.26), note that am is the number of acoustic 
modes considered in the system. The nondimensional panel shape 
function is defined as 




Note that the quantity 



is nondimensional; however, the actual form of k x ij cannot be 
determined until an appropriate choice for \|/ij has been 
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determined; w r will then be used to nondimensionalize the 
frequency appearing in the k x ij • 


.Wtinn 4.3; Recove r ing Some Familiar Rffillltfi From J hS. 
General Re *l Panel Analysis 


A number of familiar example problems may be recovered 
from the real panel governing equations, either in dimensional 
form, such as Equation (4.2.22) with appropriate substitution for 
the elements M «P pq , Jij, T«P ijpq , Q E a Ppq, and Tab ijcdf or in 
nondimensional form, such as Equations (4.2.24) through (4.2.27). 

The most basic example which can be recovered is the case 
of a single panel in an infinite duct. After appropriate substitution. 
Equation (4.2.22), the general dimensional governing equation, 
reduces to (neglecting all summation limits temporarily) 



a P P q 


-XXXXA 

c d r s 


C r1Xxf 

i j J 


4 >pq Vij dS p 


aP 


lpor 


J ^pq Vij d Sp d 

h 


dStotal 


■i 


PeI ‘t'pq dS“ p 


(4.3.1) 


For a single panel in a duct, all acoustic modes except the one 
dimensional mode may be ignored; as such, yij -> 1 a nd ^xij t°/ c * 
For a flat panel, the mode shape * -> 1; all summations can be 
removed. From Equation (4.3.1), the ratio of the integral over 
dS p cd to the integral over dS to tal -> 10 Equation (4.3.1) reduces to 

A M [o)| + 2iCa><o n - to 2 ] - iAWHpcoc = P E WH (4.3.2) 
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Note that icoA = U, the panel velocity; also, p t = pcU. M, the mass 
per unit area, may be replaced by m, the panel mass. From a fixed 
wall approximation, P E = 2Pi. These manipulations yield 

u n jtOn + 2it>o>n - to 2 ] _ _ 2 pj (4.3.3) 

WH io) 

The panel's mechanical impedance may be expressed in the 
familiar form 

Z™ = R + ijmco -^] (4.3.4) 


where R is the mechanical damping and s is the spring constant. 
Also, 

r = R®n (4.3.5) 

s 2s 

Substitution of Equations (4.3.4) and (4.3.5) into (4.3.3) yields a 
familiar relationship for the ratio of transmitted to incident 
pressure in a duct as 


2 (4.3.6) 

Pi Zm _ j 
WHpc 

A second limit that can be recovered is the branch analysis 
relationship for the ART 4-panel geometry in a duct. 27 The branch 
analysis result gives the ratio of the transmitted pressure to the 
incident pressure across the panel barrier in a duct when only the 
one dimensional acoustic mode is considered. In dimensional form 
for an anechoic termination, 

Pt = A 

Pi A + — 2 — Zn Zn Z21 Z22 

pcWH 
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(4.3.7) 



where 


A = Zn Z12 Z21 + Zi 1 Z12 Z22 + Zn Z21 Z22 + Z12 Z21 Z22 

Note that all impedances in Equation (4.3.7) are mechanical 
impedances, of the form as shown in Equation (4.3.4). Again, \|/ij > 
1, kxij -> co/c, and the panel mode shape <|> -> 1. With respect to the 

dimensional governing Equation (4.3.1), summation over the panel 
indices ap and cd from 1 to 2 is now required. Following a similar 
logic as that described to derive the governing equation for a 
single panel in a duct, Equation (4.3.1) will reduce to a system of 
four nondimensional equations of the form 


U 


11 


M Zm„ + 1 - 

ns 2 / 

4 j 

+ u, 2 

N i™, 

. \s 2 l 


Mu 

4 


Mh + LLu + u 2 i 

44 " 


+ Ui2 + U2L + M21 + P t = 2 
4 4 4 


+ M 21 + M 22 . + p T = 2 

4 4 


.B Zm2 ' + 4 


4. M 22 . + p x = 2 

4 


Un + u u + y n + u 2 2 

4 4 4 




+ Px = 2 

(4.3.8(a-d)) 


where nondimensional velocities are denoted by 


Uij 


Uij P £ 
Pi 


and nondimensional impedances are given by 
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Zjnjj — MRjj coy 




(Oij 


(D 


Additionally, mass ratios are defined as 


MRjj 




mn 


and 



Other relationships have been defined previously. Finally, a fifth 
equation expresses the fact that the resultant nondimensionalized 
pressure ratio is an average of all nondimensional panel velocities, 
that is, 


Ujj-4 ]hi +J^hi +lhi _p t = o 

4 4 4 4 

(4.3.9e) 


With the aid of the symbolic manipulation program Mathematica, 
Equations (4.3.9) can be solved for the unknowns Ull, Ul2» U21, 
U22, and PT- In nondimensional form, PT is the equivalent of the 
dimensional branch analysis relationship shown in (4.3.7) 

p t = A. — — (4.3.10) 

A + N Zn Zi2 Z21 Zm 


A = Zn Z12 Z21 + Z11 Z12 Z22 + Zn Z21 Z22 + Z12 Z21 Z22 
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Section 4.4; Structural Branch Analysis 


The general solution of the governing Equations (4.3.1) for 
the real panel analysis would formulate as described below. The 
solution would be carefully formulated for a general set of mode 
shapes; this linear system (of the generic form [A] x = [B]) would be 
sized on the order of the (number of panels) * (number of modes 
per panel). Numerical solution would be required to obtain the 
normalized modal amplitudes x, which in turn would determine 
the normalized panel modal velocities. The [B] matrix contains 
expressions composed of the exterior pressure loadings integrated 
over the particular panel limits (where the loading is applied) with 
its associated mode shape. 

None of the latter issues presents any real problem; the 
major difficulty in the exact solution to this problem is the terms 
resulting from the fluid loading of the higher structural mode 
shapes on other panel modes. The [A] matrix for the system would 
be constructed in the same manner as the governing equations for 
the Chapter 2 model were formulated; that is, the diagonal terms 
of the matrix would contain the panel modal impedances, and all 
terms would contain the generalized forces. The latter forces 
result from modal apparent mass loadings applied to each panel 
modal impedance (including the panel itself) from all structural 
mode apparent mass contributions of each plate in the array. For 
the Chapter 2 models with panels represented by a rigid flat plate, 
this procedure was somewhat tedious but nonetheless possible. In 
those problems, the flat rigid panel shapes presented an obvious 
symmetry with respect to the apparent mass loadings (generalized 
forces) which could be exploited. As shown in Figure 2.2.1 for flat 
rigid panels, the apparent mass loading on panel 12 resulting from 
the motion of panel 1 1 is the same as the loading on panel 1 1 due 
to the motion of panel 12. Such symmetry does not exist in the 
general real panel problem. Figure 4.4.1 displays this concept 
graphically, where the top diagram shows two adjacent panels 
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panel 1 


panel 2 



panel 1 panel 2 


Figure 4.4.1: Dissymmetry of apparent mass modes with 

the inclusion of higher structural modes in the real panel 
analysis. 

with pinned boundary conditions on all edges; panel 1 on the left 
has mode 22 only, while panel 2 on the right has mode 11 only. 
Now consider the generalized force on panel 2 mode 11 resulting 
from panel 1 mode 22. Looking left from panel 2, panel 1 (with 
mode 22) does not provide the same fluid loading as the situation 
shown in the lower portion of Figure 4.4.1, where the panel mode 
shapes are reversed. In the upper case shown in Figure 4.4.1, 
where upward panel deflections existed in the panel 1 mode 22 
shape on the left, the lower portion of the figure shows that now 
downward panel deflections are present instead. Rephrased, panel 
1 mode 22 does not load panel 2 mode 11 in the same manner as 
panel 1 mode 11 loads panel 2 mode 22. Thus symmetry cannot 
be exploited in the calculation of the general apparent mass 
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flat panel mode shape 



clamped panel mode shape 


Figure 4.4.2: Qualitative apparent mass differences for 

flat, pinned, and clamped mode shapes. 

loadings or generalized forces on the panels, and this is the major 
impediment to a full solution of the problem. 

Another issue worthy of discussion is the general magnitude 
of the apparent masses for the various mode shapes which might 
be used. Figure 4.4.2 shows a side view of two panels with three 
possible mode shapes. The dotted line in each sketch represents 
the panel zero displacement line, and the panel mode shapes 
sketched represent the maximum deflected panel mode shapes 
which the surrounding air mass will load. The flat mode shapes 
shown in the top portion of Figure 4.4.2 have an idealized square 
displacement envelope, where much fluid is sloshed back and forth 
because the portions of the panels butted together in the middle 
are discontinuous and separate. The pinned mode shapes are 
piecewise continuous, as shown in the middle of Figure 4.4.2, and 
this overall panel displacement envelope results in less fluid 
sloshing back and forth between the panels. Finally, the lower 
portion of Figure 4.4.2 shows adjacent clamped panel sections, 
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where the panel displacement envelope is now totally continuous; 
even less mass is sloshed back and forth between the panels. The 
model presented in Chapter 2 used rigid flat plate mode shapes, 
which represents a conservative estimate of the apparent mass 
loading effects as described graphically in Figure 4.4.2. The 
Chapter 2 model indicated that while there is some effect due to 
the panel apparent mass loading, it is not severe unless the ratio of 
surrounding fluid mass to the panel mass is high. This latter 
parameter is S 2 as shown in Tables 2.2.1 and 3.3.1; if S 2 is less than 
0.1, the apparent mass effects are not severe, and it can be 
predicted that the frequencies of maximum ART noise reduction 
will shift to only a slightly lower value due to the additional mass 
of the surrounding fluid in the vibrating panel-fluid system. 

The branch analysis method is a reduced solution method 
which can be used to obtain a model which will predict the general 
acoustic behavior of real panels in an ART configuration. This 
approach will use only the one dimensional acoustic mode and as 
such, ignores apparent mass loading. However, it will permit an 
investigation of the effects of higher structural modes on ART 
tuning. The analysis follows a schematic identical to Figure 2.2.2, 
except now more branches are required to model the individual 
structural modes considered on each panel. A noise reduction in 
decibels may be described as 


NR = -10 log 



2 

-10 log 


P B , 

LMJ 

P B 


(4.4.1) 


where 


P B Zb cos k(L - Lm] + ip 0 c sin k[L - Lm] (4.4.2) 

d b Zb cos kL + ip 0 c sin kL 

v 0 
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Note that these quantities have been defined in Chapter 2 and 
carry similar meaning here. Note that Equation (4.4.2) describes 
the duct acoustics exclusive of the panel wall barrier (also 
equivalent to the branch analysis of Chapter 2). The remaining 
analysis portion requires the general expression of the pressure 
ratio P 0 B /Po A to include higher structural modes. Using the branch 
analysis procedure explained in Chapter 2, that pressure ratio may 
be determined as 


where 


P B 

iJL = 
P A 

Pq 


2 2 modes modes 

2 2 2 2 

a=l P=1 p=l q=l 


ia j 


4$ dS«|> 


c 2 2 modes modes 

x x 

zj a=1 P =1 P =1 



(4.4.3) 



Zb cos kL + ip 0 c sin kL 
p 0 c cos kL + iZb sin kL 


(4.4.4) 


Sa is the overall cross sectional area of the duct, dS a p corresponds 
to integration over the area defined by panel ocp. A generalized 
area coefficient may be described as 



(4.4.5) 


where <J> a P pq is the panel mode shape for panel ap, mode pq, S a p is 
the area of branch ap, and 
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M “5 = < q [(a“P * '«)L, + 2tS < (a“5 * + <’ K e '«} 


(4.4.6) 


where the usual convention follows that a and P refer to the 
matrix-like panel identification number (see Figure 2.2.1) and p 
and q index the individual panel modes. Note that a, p, p, and q 
again apply to the mass, natural frequency, modal amplitudes, and 
damping ratios of Equation (4.4.6) The generalized mass M a Ppq in 
Equation (4.4.6) above may be defined as 28 


< = < 



SaP 


(4.4.7) 


Note that nondimensionalization removes the modal amplitude 
coefficients shown in Equation (4.4.6) by converting them to 
nondimensional velocities as 


icoA^q = U^q (4.4.8) 

Using the nondimensional parameters defined earlier in this 
chapter and in Table 2.2.1, the results above may be expressed in 
nondimensional form as 


P B 

fJL. — 

pA 

^0 


2 2 modes modes 

III I 

a=l P=1 p=l q=l 


icoS2 j 


— ap — 
<|>pq dSap 


4SiA/pq 


, 2 2 modes modes 

J- + l l I I 

a=l |)»1 P=1 <1=1 



4SiMpq 


(4.4.9) 
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where integrations are now calculated over nondimensional panel 
limits between -1 and +1. Also, 


-«P 

pq 


A/; r n = M 


a p 

pq 


-ap 2 



(4.4.10) 


g Z h cos ({ 0 S 1 S 4 ) + i sin (( 0 S 1 S 4 ) 
cos (C 0 S 1 S 4 ) + iZb sin ( 00 S 1 S 4 ) 


(4.4.11) 


and 


P® L= Z h cos(toS 1 (S 4 - S 3 )) + i sin(coSi(S 4 - S 3 )) (4 4 12 ) 

pB Zbcos (©S 1 S 4 ) +i sin (©S 1 S 4 ) 

Finally, the choice of panel mode shapes must be made. For flat 
panel mode shapes, _ 

Pinned panel mode shapes are represented by 


<Xy,z) = sin (prcy) * sin (qjtz) 


with higher mode nondimensional resonance frequencies for 
simply supported square panels analytically denoted by ^8 

5„ = miAj£ 

2 

The presence of higher structural mode shapes for the pinned 
cases are added to the model in the order mode 11 , mode 12 , mode 

22 , mode 13, mode 23, and mode 33. 

Clamped panel mode shapes may be represented according 

to Young by 30 


<Ky»z) = 


p q _ _ 

S X ^mn Xm (y)Xn (z) 
m=l n=l 
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where 


X(y) = [cosh(e m y) - cos (e m y) - a m (sinh (e m y) - sin (e m 


X(z) = [cosh(e n z) - cos (e n z) - a n (sinh (e n z) - 



Natural frequencies, amplitude coefficients, £m* an d schematic 
panel mode shapes are as given in Figure 4.4.3 (reproduced from 
Leissa 31 )- Depending on the choice of mode shapes, analytic or 
numerical integration is used to determine the nondimensional 
weighting coefficients defined as C a rca and CmassI these are 
denoted as 


Carea“q = 


i 


“up 

<|>pq dS<xp 


Sap 


(4.4.13) 


Cmass“q = 


.[[ 


-nap 

0pq dSap 


Sap 


(4.4.14) 


It is important to note that this analysis says nothing about 
how many modes will be present on a given panel; that latter 
quantity is judiciously chosen when performing noise reduction 
calculations. For small panels such as those used in the analyses in 
Chapter 2 and 3 (aluminum panels 12" by 12"), only a few modes 
(and perhaps essentially only one mode!) might be present. Much 
larger panels might have additional modes present. 
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Figure 4.4.3: Frequency parameters, amplitude coefficients, and schematic mode shape 

configurations showing modal lines for clamped modes as determined by Young (reproduced from 
Leissa). 








Section 4.5: Structural Branch Analysis ReSttllS 


Section 4.5.1 1 Mode S tudy and Co mpar i son of ART a M 

identical — Panels 

As shown in Figure 4.4.3, Young determined mode shapes 
and natural frequencies for the first six modes of a clamped plate. 
While it is probably not feasible to expect six modes to exist on a 
square plate 12" by 12", the simulations were performed for all six 
modes for both clamped and pinned geometries to ascertain 
possible performance. Damping ratios C for a11 modes of a11 P lates 
are set to a value of 0.01. Figure 4.5.1 through 4.5.6 show noise 
reduction predictions for these 6 cases. Figure 4.5.1 makes a 
comparison of identical and ART tuned panel noise reductions for 
flat (rigid) panel shapes (as predicted by the Chapter 2 branch 
analysis), pinned, and clamped mode shapes. The general trend 
here is that the greatest noise reduction is observed with clamped 
plates (presumably because they are stiffer); flat plates provide 
the least amount of noise reduction. The noise reduction 
predictions generally show about 3 dB noise reduction increase for 
pinned panels over flat rigid panels, and again about 3 dB noise 
reduction increase for clamped plates over pinned plates. These 
differences are observed at both low and high frequency limits. 
Thus the clamped panel appears more mass and stiffness 
dominated than the pinned plate, and the pinned plate is equally 
varied in stiffness and mass domination from the flat rigid plate. 
Another noticeable difference is that the ART frequency 
bandwidth at a constant noise reduction widens from flat to 
pinned to clamped plates. This widening is a typical result seen 
for reduced apparent mass loadings with the two panel and four 
panel Chapter 2 models. This noise reduction bandwidth increase 
results because less apparent mass is sloshed back and forth 
between the panels since less panel surface is moving. Fixed 
adjacent edges of the real panels also contribute to less apparent 
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Figure 4.5.1: Noise reduction prediction for ART and 

identical panels using 1 structural mode and panel mode 
shapes as indicated. 



Figure 4.5.2: Noise reduction prediction for ART and 

identical panels using 2 structural modes and panel mode 
shapes as indicated. 
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mass sloshing as depicted graphically in Figure 4.4.2. 

Consequently it can be expected that the actual ART cancellation 
frequency will not be decreased (apparent mass loaded to reduce 
the ART noise reduction maximum) as much as indicated by the 
Chapter 2 model. 

Figures 4.5.2 and 4.5.3 show noise reduction predictions for 
pinned and clamped panels with two and three structural modes 
acting. The value of the weighting coefficient C a rea described by 
Equation (4.4.13) is zero for both the second and third structural 
modes for both clamped and pinned mode shapes, and as such 
these modes play no role in changing the noise reduction 
predictions. It is only when the fourth mode (mode 13) is added 
to the pinned calculation that any change is seen, and this change 
occurs in the identical panel calculation in Figure 4.5.4. Mode 13 
has nonzero weighting coefficients, and the natural frequency is 
five times the mode 11 natural frequency. Thus a single plate 
section has effectively two resonance frequencies, and the ART 
effect can occur. This mode 13 natural frequency occurs at co = 5.0, 
and an ART cancellation is seen at w ~4.75. In the ART panels case, 
the lowest ART panel natural frequency is 0.5, and five times this 
frequency yields co = 2.5, which is another panel resonance 
frequency, so no net noise reduction effect is observed. The fourth 
mode for clamped plates also has a zero area weighting coefficient; 
no additional noise reduction contributions occur here. With the 
addition of the fifth mode (mode 23 for pinned and mode 5 as 
shown in Figure 4.4.3 for clamped), clamped panels also have 
higher modes playing an ART role, and additional noise reduction 
peaks are seen in Figure 4.5.5. The sixth structural mode (mode 
33 pinned and mode 6 clamped as per Figure 4.4.3) plays no 
significant effect. Therefore, while it is theoretically possible to 
show the positive ART influence of higher structural modes for 
both pinned and clamped panel shapes, these ART contributions 
would be occurring at very high frequencies. Higher mode effects 
would be quite minimal, if existent at all. 
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Nondimensional Frequency 


Figure 4.5.3: Noise reduction prediction for ART and 

identical panels using 3 structural modes and panel mode 
shapes as indicated. 



Nondimensional Frequency 

Figure 4.5.4: Noise reduction prediction for ART and 

identical panels using 4 structural modes and panel mode 
shapes as indicated. 
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Figure 4.5.5: Noise reduction prediction for ART and 

identical panels using 5 structural modes and panel mode 
shapes as indicated. 



0 1 2 3 4 5 

Nondimensional Frequency 


Figure 4.5.6: Noise reduction prediction for ART and 

identical panels using 6 structural modes and panel mode 
shapes as indicated. 
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Figure 4.5.7: ANNR and ART noise reduction predictions 

using two structural modes and panel mode shapes as 
shown . 


Finally, Figure 4.5.7 shows a comparison of noise reduction 
predictions using ART tuned pinned and clamped panel sections 
compared to the ANNR calculation for identical pinned and 
clamped panels. This prediction is very similar to previous ART 
and ANNR comparisons; clearly ART tuned panels are performing 
much better with at least a 15 to 20 dB noise reduction increase. 

Section 4.5 .2: Parameter Shldfcs 

Figures 4.5.8 through 4.5.11 show relative parameter study 
results; these results are very similar to those obtained using the 
Chapter 2 model. Figures 4.5.8a and 4.5.8b show a parameter 
study with increasing apparent mass loading for pinned and 
clamped mode shape panels, respectively. As reported earlier, 
increasing apparent mass loading by successive doubling of this 
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Nondimensional Frequency 


Figure 4.5.8a: Parameter studies varying the panel 

apparent mass loading as indicated for 1 structural 
mode with a pinned mode shape. 



Nondimensional Frequency 

Figure 4.5.8b: Parameter studies varying the panel 

apparent mass loading as indicated for 1 structural 
mode with a clamped mode shape. 
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Figure 4.5.9a: Parameter studies varying the panel 

damping as indicated for 1 structural mode with 
a pinned mode shape. 



0 12 3 4 

Nondimensional Frequency 


Figure 4.5.9b: Parameter studies varying the panel 

damping as indicated for 1 structural mode with 
a clamped mode shape. 
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Nondimensional Frequency 

Figure 4.5.10a: Parameter studies varying the ART 

noise reduction bandwidth by increasing the higher 
panel natural frequency for 1 structural mode with 
a pinned mode shape. 



Nondimensional Frequency 

Figure 4.5.10b: Parameter studies varying the ART 

noise reduction bandwidth by increasing the higher 
panel natural frequency for 1 structural mode with 
a clamped mode shape. 
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Nondimensional Frequency 

Figure 4.5.11a: Parameter studies varying duct termination 

impedance as indicated for 1 structural mode with a pinned 
mode shape. 



0 1 2 3 4 5 


Nondimensional Frequency 

Figure 4.5.11b: Parameter studies varying duct termination 

impedance as indicated for 1 structural mode with a clamped 
mode shape. 
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parameter value tends to narrow the ART noise reduction 
bandwidth in more or less uniform amounts. 

Panel damping ratios are varied as shown in Figures 4.5.9a 
and 4.5.9b for pinned and clamped panels, respectively. In both 
pinned and clamped cases, as the damping ratio is increased, noise 
reductions at panel natural frequencies are increased due to less 
acoustic transparency. Maximum noise reductions at the ART 
design frequencies are reduced as shown. However, the noise 
reduction bandwidth at about 25 dB has remained essentially 
constant for all damping ratios. This prediction implies that the 
ART concept will behave in a robust manner when applied with 
real panels sections. 

Figures 4.5.10a and 4.5.10b show increasing ART noise 
reduction bandwidth by increasing the higher natural frequency 
panel resonance frequency as shown for pinned and clamped 
panels, respectively. Note also that the maximum amount of noise 
reduction obtained has increased as well. Increasing the higher 
natural frequency panel resonance frequency always provides 
better ART performance. 

Figures 4.5.11a and 4.5.11b show variations in the duct 
termination impedance as indicated. As shown in Chapter 2, any 
termination impedance harder than Zb = 10 + Oi produces about 
the same noise reduction prediction. “Hardbox” frequencies again 
occur as predicted by Equation (2.4.4), and the ART method still 
provides noise reduction even with the presence of a reflected 
interior pressure wave from a duct termination impedance. 
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Chapter § 

Experimental C onfirmation of thfi — ARX — CflllCCtt i 


section 5.1l introdiKtim 

A series of experiments has confirmed the potential of 
Alternate Resonance Tuning for noise reduction applications in 
simple paneled structures. This chapter presents highlights of the 
various experiments conducted. 

Section 5.2i Experimental — Salim 

Figure 5.2.1 shows an overall schematic of the experiment. A 
sound source driven by a pure tone generator or a white noise 
generator is coupled to a duct. A panel array test section is placed 
between the source-duct coupling and the downstream portion of 
the duct. Microphones monitor the overall sound pressure levels 
in decibels (dB) upstream of the panel test section (closer to the 
sound source) and downstream of the panel test section as shown 
in Figure 5.2.1. Additionally, a microphone can be placed in the 
center of the panel test section if the test section is a double wall 
configuration. Termination impedances available include anechoic 
and hard wall impedance values. Data shown herein considers 
only the anechoic termination impedance, since the ART effect is 
most clearly observed without the interference of a reflected wave 
from the downstream termination. 
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Figure 5.2.1: Overall experimental setup showing general 

microphone positions, test section location, and processing 
equipment. 
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Two signal processing methods are used. One setup requires 
the use of a white noise generator and a Scientific Atlanta SD-380 
spectrum analyzer; in this mode of operation, both microphone 
signals are fed to the SD-380 operating in the transfer function 
mode to produce the noise reduction across the panel barrier. The 
second procedure is a modern variation on a spectrum analysis 
method most often used by the Bruel & Kjaer company, where a 
pure tone generator is mechanically driven through the frequency 
range of interest. Frequency is monitored directly from the pure 
tone generator. The microphone signals are filtered and converted 
to a direct current signal which is linearly proportional to sound 
pressure level. A Macintosh II computer then processes the 
signals using LabVIEW, an icon-based programming tool. 

Lab VIEW collects the DC voltage signals in integer data form and 
processes them using appropriate calibration to generate 
frequency and sound pressure level information. A noise 
reduction is then calculated by differencing the upstream and 
downstream microphone levels. The LabVIEW software creates 
almost immediate graphical representation of the results; 
additionally, the program can write frequency information and 
individual microphone levels to a data file for later processing. 

Panels used in the panel arrays are (with one notable 
exception mentioned later) passively driven audio speakers. These 
“panels” are low in cost and have resonance frequencies which can 
be easily modified with the addition of mass to the cone, lowering 
natural frequencies. Epoxy can also be carefully added to the 
suspension of the panel to raise the stiffness and increase the 
natural frequency. In this manner, the required separation of the 
natural frequencies of panels of equal size can be obtained, and the 
panels can be forced to work against each other acoustically 
according to the ART principle. The voice coil on the speakers 
provides a handy transducer for measuring the phase difference 
between the panels. 
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As mentioned earlier, there is one exception to the use of 
passively driven audio speakers as panel sections; Section 5.7 
presents data obtained using real panel sections composed of 
Bakelite. 

Section 5.3i single Panel Wall Noise Reductio n 
Measurements 

The speakers generally used for the single panel wall 
experiments are 4" in diameter and have an unmodified natural 
frequency of about 200 Hertz. These “panels” are arranged in the 
duct in a number of ways as shown in Figure 5.3.1, with single 
panel, double panel, and four panel wall configurations possible. 
After modification, ART panels have resonances of approximately 
100 Hertz and 300 Hertz; use of other natural frequencies will be 
noted as required. 

Figure 5.3.2 shows noise reduction measurements for 
identical panels and ART tuned panels. Both panel types exhibit 
approximately stiffness-dominated behavior at low frequencies 
(well below 100 Hertz) and mass dominated behavior at high 
frequencies (above 400 Hertz). The identical panel results show a 
noise reduction minimum at the panel natural frequency of 200 
Hertz. The uneven noise reduction decrease below the natural 
frequency and uneven noise reduction increase above the natural 
frequency is due to the lack of a perfectly anechoic termination. 

(A pure anechoic situation is quite difficult to simulate in a 
laboratory.) The ART tuned panels show minima at the natural 
frequencies of 100 and 300 Hertz; between these frequencies, a 
maximum noise reduction of about 24 dB is observed. The 
proceeding results indicate the robustness of the ART method. 

The phase difference between the panels can be measured 
using the speaker voice coils, as shown in Figure 5.3.3. Panels are 
in phase at low frequencies up to the lowest panel natural 
frequency. Between the panel natural frequencies, a phase 
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Figure 5.3.1: Single wall panel array test sections. 
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Figure 5.3.2: Experimental noise reduction measurements 

with two ART tuned panels and two identical panels. 
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Figure 5.3.3: Measured phase difference between two 

ART panels used to generate ART data shown in 



Figure 5.3.2. 
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difference between panels of about 165 degrees was noted. At 
high frequencies, the panels are again in phase. Figures 5.3.2 and 

5.3.3 demonstrate the possibility of tuning panels to achieve the 
ART requirements of out-of-phase behavior and approximately 
equal amplitude motions. 

The concept of ANNR can be applied to experimental data to 
determine the true ART contribution to noise reduction. Figure 

5.3.4 shows the two panel ART data and an ANNR calculation 
obtained from individual panel noise reduction data*. The ART 
concept effectively adds over 10 dB of noise reduction when 
compared to the averaged noninteracting noise reduction 

calculation for individual panels. 

Figure 5.3.5 shows experimental noise reduction data with 
the two panel setup obtained by varying the higher panel natural 
frequency as shown. Note that increasing the frequency difference 
between the high natural frequency panel and the low natural 
frequency panel increases the noise reduction bandwidth as well 
as the maximum noise reduction value obtained. In general, this 
strategy proves advantageous from both theoretical and 
experimental considerations. 

Four panel tuning is shown experimentally in Figure 5.3.6, 
where ART panel natural frequencies are arranged at about 90, 
190, 280, and 410 Hertz. While these are certainly not optimum 
tuning frequencies, the data does indicate the possibility of tuning 

for more than one ART design frequency. 

It is possible to verify experimentally the exponentially 
decaying nature of the sound field in the downstream vicinity of 
the panel wall using a modified impedance tube setup as shown in 
Figure 5.3.7. A rolling cart microphone with a long extension tube 
is placed through the sound source, duct coupling, and panel wall 


* ANNR data is obtained by measuring noise reduction using single panels m 
the single panel duct. A separate noise reduction data set is then obtained 
for the 100 Hertz and the 300 Hertz natural frequency panels, and Equation 
(2.4.1) is invoked to determine the ANNR. 
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Figure 5.3.4: Two panel ART noise reduction measurement 

plotted with two panel ANNR calculation. 



as shown. 
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400 Hertz. 


at a position as shown in the panel array cross section detail in 
Figure 5.3.7. Note that the end of the microphone tube is 
supported on a rolling wheel, insuring measurement at the same 
vertical position in the duct. This apparatus allows the 
measurement of sound pressure level falloff at a given frequency 
at various longitudinal positions in the duct relative to the 
coordinate system as shown. Figure 5.3.8 displays the sound 
pressure level difference in decibels between the upstream 
microphone and the cart tube microphone plotted as a function of 
distance down the duct at various frequencies. Note that at the 
100 Hertz frequency (lower panel resonance frequency) there is 
little sound level falloff; recall that at the natural frequency, a 
panel will appear acoustically transparent in the limit of no 
damping. However, at 200 Hertz (the frequency of maximum noise 
reduction), an 18 dB sound pressure level drop has occurred across 
the panel wall; an additional 6 dB of noise reduction occurs by 
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FALLOFF (dB) NOISE REDUCTION (dB) 



Figure 5.3.8: Difference between upstream microphone and cart 

microphone levels as a function of distance down the duct for the 
two panel ART array. Measurement is taken as shown in Figure 
5.3.7. 



DISTANCE (centimeters) 

Figure 5.3.9: Falloff measured away from the panel wall array 

with the cart microphone. The falloff in decibels is referenced 
to the sound pressure level just downstream of the panel wall. 
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moving approximately 25 centimeters away from the panel wall. 
Sound pressure levels near the wall are higher; however, ART 
tuning has reduced them significantly. Figure 5.3.9 shows sound 
pressure level relative to the level just downstream of the panel 
array; at 200 Hertz, the exponential decay is much more evident 
on this scale. 

Section 5.4: Double Panel Wall Noise Reducti on 
Measurements Wi th One Identical Wall and One ART Wal l 

Double wall structures are easily investigated with the test 
section shown in Figure 5.4.1. An upstream panel array holds two 
identical panels with resonance frequencies of 200 Hertz; in the 
downstream wall, identical or ART tuned panels may be arranged. 
Figure 5.4.2 shows experimental noise reduction results for 
identical downstream panels (meaning all four panels have 
resonance frequencies of 200 Hertz) and ART tuned downstream 
panels. Since the panel test section now contains two panels in 
series, the noise reduction falloff at low frequencies is about 12 
dB/octave in the low frequency limit (stiffness dominated 
behavior). Noise reduction increase in the high frequency limit 
should theoretically be 12 dB/octave (mass dominated behavior). 
The low frequency and high frequency noise reduction results in 
Figure 5.4.2 show the effect of these higher order impedances. The 
noise reduction at the panel natural frequencies also increases 
because the sound must now travel through two panel arrays. 

ART tuning has achieved about 18 dB of additional noise reduction 
at 200 Hertz. Figure 5.4.3 shows the measured phase difference 
between the downstream panels for this same case; higher order 
impedances make this phase difference data appear more 
“rounded” than the phase difference data for the single array 
(Figure 5.5.3). 

The effect of panel wall spacing is shown in Figure 5.4.4 for 
panel wall array coupling lengths of 7", 2", and 1". High frequency 
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Figure 5.4.1: Top view of double wall test section with identical 

panels upstream and ART tuned or identical panels downstream. 
Midstream microphone located as shown. 
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Figure 5.4.2: Experimental noise reduction measurements 

in the double wall setup with four identical panels and 
identical panels upstream and ART tuned panels downstream 
with the 7" coupling. 



Figure 5.4.3: Phase difference between ART tuned 

panels used in noise reduction data shown in 
Figure 5.4.2. 
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Figure 5.4.4: Noise reduction measurements for ART tuned 

panels downstream and identical panels upstream with 
coupling separation distances as shown. 


data for the 7" spacing shows more noise reduction than the 
narrower spacings since the larger spacing has permitted the two 
panel arrays to behave more independently. In the ART 
frequency regime, the noise reduction data are not drastically 
different either in noise reduction bandwidth or maximum noise 
reduction obtained. 

Wtinn 5.5: Double Pa nel Wall NoiSff Reduction 
Measurements With Tw ft ART Walls 

Double wall geometries permit the investigation of noise 
reduction results for two sets of ART tuned panels, and Figure 
5.5.1 displays the arrangement possibilities. Inline arrangement 
places panels of identical natural frequency on the same side of 
the duct; opposed arrangement places panels with the same 
natural frequencies on opposite sides of the panel array. The 
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Figure 5.5.1: Double wall test sections with ART panels upstream 

and downstream. Note difference between inline and opposed 
ART configurations. Midstream microphone located as shown in 

Figure 5.4.1. 

major physical difference here is the loading of the higher 
evanescent modes or surrounding fluid; inline panels slosh fluid 
back and forth lengthwise in the duct, whereas in the opposed 
geometry, panels tend to slosh fluid in a crosswise manner. Figure 
5.5.2 indicates that the results are quite similar for a 7 coupling 
spacing except at the maximum noise reduction, where the 
opposed panel geometry appears about 8-10 dB better. 

Examination of noise reduction results in Figure 5.5.3 for inline 
panel arrangements at various spacings show that relatively close 
spacings of 1" or less behave in a similar manner; the 7' spacing 
shows increased noise reduction behavior. All noise reduction 
maxima occur at about the same frequency in Figure 5.5.3, 
indicating that the fluid loading on the inline panels is at least 
approximately constant. However, in the opposed ART panel 
results shown in Figure 5.5.4, successively closer spacings lower 
the frequency of maximum noise reduction, indicating that the 
panels may have an effectively higher fluid loading at closer 
spacings. 
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Frequency (Hertz) 

Figure 5.5.2: Experimental noise reduction measurements 

for inline and opposed double ART walls. 



Frequency (Hertz) 

Figure 5.5.3: Experimental noise reduction results for 

inline double ART walls with coupler spacings as shown. 
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Figure 5.5.4: Experimental noise reduction results for 

opposed double ART walls with coupler spacings as shown. 


Section 5.6: Double Panel Wall Noise Reductio n 
Measurements Wi th Noniiniform Panel — Si££ 

Figure 5.6.1 displays another panel arrangement possibility 
tested in the four panel array setup, where the double wall 
consists of one large upstream panel (free standing resonance 
frequency of 150 Hertz and 8" diameter) and four small 
downstream panels. Noise reduction data for identical 
downstream panels and ART tuned downstream panels are 
displayed in Figure 5.6.2, where ART tuning results in a noise 
reduction increase of about 22 dB over identical tunings. These 
results show the effect of the larger panel with a lower resonance 
frequency; for identical panels downstream, the large panel 
produces a noise reduction minimum at about 115 Hertz, for ART 
tuned panels, this minimum occurs at about 130 Hertz. ART panels 
apparently reduce the fluid loading on the larger panel; as a result, 
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Figure 5.6.1: Double wall test section with nonuniform panel sizes. 

ART or identical panels are placed on downstream panel. Midstream 
microphone is placed on duct centerline and equidistant from each 
panel array. 
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Figure 5.6.2: Experimental noise reduction comparison 

between uniform and ART panel arrangements for the 
1" coupling. 



Figure 5.6.3: Noise reduction data from double wall 

nonuniform panel size experiments with coupling 
spacings as shown. 


173 



the frequency of minimum noise reduction associated with the 
large panel natural frequency increases from the identical to the 
ART panels case. Figure 5.6.3 compares ART tuning data for three 
coupling spacings; for the largest spacing shown (2"), the ART noise 
reduction peak is the narrowest; all coupling spacing data shows 
very similar behavior in and around the ART frequency regime. 
Relative maxima and minima at frequencies above 300 Hertz show 
the effect of coupling resonance, and the largest separation data 
shows a more distinct attempt to achieve mass dominated, high 
frequency behavior at a lower frequency. 

Section 5.7: Noise Reduction MeaSUrfmfTltS With Re al 
Panel ART Walls 

Experimental results presented to this point have been 
obtained using passively driven speaker assemblies as panel 
sections. The assembly shown in Figure 5.7.1 allows the 
experimental investigation of the ART effect using real plates. 
Multiple socket head screws are used to clamp Bakelite plates 
(1/64" thickness) into an aluminum frame for use in the 
experimental apparatus shown in Figure 5.2.1. Figure 5.7.2 shows 
the noise reduction results for a plain panel with no modifications 
with a solid line. This data shows stiffness dominated, 6 dB noise 
reduction rolloff behavior at low frequencies. Again the 
irregularities in the low frequency data indicate the lack of a 
perfect anechoic termination. The panel natural frequency occurs 
at about 475 Hertz. Other small peaks in the noise reduction data 
are due to spurious resonances in the duct at higher frequencies. 

It is useful to compare the experimental data with plate 
theory to determine just how predictable the panel behavior is. 
Panel noise reduction behavior has been observed in the 
laboratory to be highly dependent on the boundary conditions. For 
example, the data shown in Figure 5.7.2 cannot be easily 


174 



Figure 5.7.1: Plan view of real panel wall test configuration. 

Noise reduction measurements are made using both single 
panel and double panel ducts. 

duplicated if the panel is not firmly attached and evenly seated in 
the frame. Any foreign matter interfering in the panel/frame 
interface can cause anomalies in the data. For the case of Figure 
5.7.2, however, the comparison between theory and data is quite 
good. For a square plate with four clamped boundary conditions, 
the natural frequencies are given by 32 

JIT (5.7.1) 

' u2 V m 


where b is the length of a side of the plate. D is the flexural 
rigidity, denoted by 


(5.7.2) 








Frequency (Hertz) 

Figure 5.7.2: Experimental noise reduction measurement 

comparison for a single 1/64" unaltered Bakelite panel 
and mass-altered Bakelite panel. 



Frequency (Hertz) 

Figure 5.7.3: Experimental noise reduction measurement 

with two panel ART setup with unaltered 1/64" Bakelite 
panel and 1/64” mass-altered Bakelite panel. 
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where E is Young’s modulus, h is the plate thickness, m is the mass 
per unit area of the plate, and v is Poisson's ratio. For Bakelite, v 
was assumed to be 0.3. are the solutions of the eigenvalue 
problem, and are shown below for the first three natural 
frequencies, along with other important values. 

~4 inches 

0.6 - 1.0 x 10 10 Pa 
0.516 kg/m 2 for 1/64" panel 

36.0 
73.8 

109.0 

Substitution of pertinent data into Equations (5.7.2) and (5.7.1) 
yields a first mode natural frequency of 496 Hertz, very close to 
the data shown in Figure 5.7.2, where a natural frequency of about 

475 Hertz was indicated. 

However, the addition of a point mass to the 1/64" Bakelite 
panel causes the panel behavior to be quite different than that 
which could be predicted by Equation (5.7.1). Figure 5.7.2 shows 
this data with a dashed line. It is curious to note the ratio of the 
two frequencies corresponding to the noise reduction minima; i.e., 

“L = &1Q=2.5 
o> 2 315 

This ratio is similar to the ratio of natural frequencies for con and 
con for a pinned plate. For a square plate with pinned boundaries 

on 4 sides 33 . 


b 

E 

m 

Xi 

*■2 

^3 



(5.7.3) 
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Plugging in m=l and n=l for con, m=l and n=2 for ( 012 , a ratio of 
natural frequencies is calculated as 


= li±2* = 2.5 (5.7.4) 

C0 H j2 + 1 2 

Arguably, the panel with additional mass may behave more like a 
pinned panel, according to plate theory equations. However, the 
presence of this second mode has created an amazing noise 
reduction between the two natural frequencies; a reduction in the 
transmitted noise of about 44 dB is seen at just below 500 Hertz. 
Recall that this is for a single panel in a single panel duct. It may 
therefore be possible to use higher structural mode behavior to 
some advantage in ART tuning a real panel. Also, the point mass 
placed at the center of the panel may help to enforce higher 
structural mode vibration. 

Figure 5.7.3 shows the noise reduction data for both the 
1/64" panel with no added mass and the 1/64" panel with added 
mass alongside each other, clamped in the two panel duct. The 
noise reduction minimum at 315 Hertz corresponds to the lowest 
natural frequency of the panel with added mass; the noise 
reduction minimum at 450 Hertz corresponds to the ©il natural 
frequency of the panel with no added mass, and an ART noise 
reduction effect of about 40 dB is observed at 350 Hertz. 

Similarly, the minima at 750 Hertz (dropped about 80 Hertz from 
Figure 5.7.2 presumably due to fluid loading) might correspond to 
the co 12 natural frequency of the panel with added mass, and 
another ART noise reduction peak of 36 dB is observed at about 
660 Hertz. If the panel with no added mass is actually behaving 
with two structural resonance frequencies, then the observed 
behavior shows two noise reduction minima as a result of three 
panel natural frequencies (two occurring from the same panel), as 
the ART theory has suggested. 
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Chapter 6 

Summary and Conclusions 


Section 6.1; Introduction 

Existing noise reduction techniques for paneled structures 
have been shown to perform reasonably well at higher 
frequencies, but are less effective at lower frequencies, 
particularly if the low frequency noise problem has high forcing 
levels such as those found in helicopter or propeller aircraft. This 
research effort has focused on a new method of noise reduction for 
paneled structures called Alternate Resonance Tuning. The 
ART concept blocks sound transmission in a frequency band by 
having adjacent panels tuned to resonance frequencies alternately 
set above and below the frequency to be cancelled. The separation 
of the panel resonance frequencies causes the panels to work 
against each other structurally and acoustically. Ideally, adjacent 
panels are designed to oscillate with equal magnitude, but are 
nearly out of phase. The problems analyzed herein have shown 
the ART concept to be of potential value in the reduction of 
transmitted noise through paneled structures. 

Section 6.2; The Four Panel Problem 

The four panel problem presented in Chapter 2 analyzed a 
periodically repeating structure of infinite extent. The building 
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block for this structure is a physically two dimensional four panel 
subsystem (three dimensional in velocity in the panel near-field) 
incorporated into a duct; the panels within this subsystem can be 
tuned to as many as four distinct resonance frequencies, damping 
ratios, and masses. An acoustic branch analysis was introduced as 
an extremely efficient tool for predicting basic behavior of this 
panel system; however, the loading of the evanescent modes was 
neglected in this analysis. A more complete analysis of normal 
incidence acoustic transmission assumed pressures and velocities 
which satisfied the reduced three dimensional wave equation. 

Using a Fourier series to model periodic displacement of a single 
rigid panel in the subsystem, an expression was derived for an 
associated velocity boundary condition which matched the panel 
displacement. This velocity boundary condition assumed that the 
one dimensional acoustic mode propagated; all other modes 
decayed due to the cutoff effect. Using the momentum equation, a 
corresponding pressure solution could be obtained; this pressure 
solution was used to determine the forces acting on the panel. 
Problem symmetry and superposition could be used to obtain 
solutions for the case where all panels were in motion. The panel 
velocities then permitted the calculation of pressures and noise 
reductions in decibels. 

Results were obtained which showed appreciable noise 
reductions for panel subsystems with two and four panels; the 
number of distinct noise reductions obtained was always one less 
than the number of distinct panel tunings within the subsystem. 
The concept of the Averaged Noninteracting Noise Reduction 
(ANNR) was introduced in order to calculate a true noise reduction 
gain due to the ART tuning. This model indicated that for typical 
system parameters, 10 to 30 dB of noise reduction could be 
obtained. Mass ratio optimization was shown to be effective in 
placing the maximum noise reduction frequencies at the desired 
values in the presence of fluid loading. A parameter study showed 
the effects of damping, panel fluid loading, ART bandwidth as a 
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function of panel natural frequencies, and variations in 
termination impedance; in all cases, the theory indicated that a 
reasonable increase in noise reduction could always be obtained by 
ART tuning. 

Section 6.3; External Pressure Field Modeling 

The purpose of this investigation was to determine the 
effectiveness of the ART concept under an external propagating 
pressure field such as that which might be associated with 
propeller passage by an aircraft fuselage. For this problem, a one 
dimensional repeating panel subsystem (two dimensional in 
velocity in the panel near-field) of N panels was analyzed, and all 
acoustic modes (both propagating and cutoff) were contained in 
the solution. The solution procedure used was similar to that 
described above for the four panel problem. A general numerical 
solution allowed the modeling of many effects, and a dual case 
study was undertaken to investigate the performance of the ART 
tuning method under typical limiting cases; these limiting cases 
included aircraft with long (low blade passage frequency) and 
short (high blade passage frequency) external propagating 
wavelengths relative to the panel size. A mode convergence study 
determined the appropriate number of modes to use for a 
converged solution. 

The results obtained indicated that ART tuning yielded noise 
reduction increases of 10 to 25 dB relative to the ANNR 
calculations. An exhaustive set of parameter studies revealed 
some very interesting effects. Avoiding mode cut on was 
determined to be an important factor in obtaining the best noise 
reduction results. Splitting the panel natural frequencies always 
proved beneficial; reasonable variations in damping ratios showed 
the robustness of the ART tuning method. Forward and reverse 
sweep direction (corresponding to different rotational directions of 
the aircraft propeller) were found to have dramatic influence on 
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the noise reduction results. This acoustic solution yielded pressure 
and noise reduction calculations covering the entire analysis 
domain, and it was possible to examine variations in a given 
parameter over the complete model area. 

A comparison was made to the model developed in Chapter 2 
to determine the effect of including all acoustic modes in the 
solution. With the exception of some deviation around mode cut 
on frequencies, it was shown that the model assumptions of 
Chapter 2 were justifiable using the Chapter 3 model. Four panel 
tuning was shown to be a distinct possibility with this geometry, 
but in the high blade passage frequency case, the effects were less 
dramatic; ART tuning was determined to be most effective at low 
frequencies. Additional miscellaneous effects such as panel 
detuning (where a realistic deviation of panel properties was 
assumed) and off-tuning considerations (where an incomplete ART 
pair was present) showed that the tuning of a fluid-structure 
interaction such as the one modeled in Chapter 3 is indeed 
complex, but that the ART concept can be used to provide 
reasonable gains in noise reduction for such a dynamic system. 

Section 6.4: The General ART Panel Analysis 

The ART panel analysis presented in Chapter 4 derives the 
governing equations for the most general case of panel geometries. 
The analysis begins by assuming general acoustic velocity and 
pressure expressions which are related by the momentum 
equation. A wall displacement function expresses the panel shape 
in a most general way; the acoustic and structural equations are 
then linked at the panel/fluid interface. Lagrange's equations are 
used to represent the equations of motion for the panels. The 
generalized forces resulting from the motion of any mode of any 
panel on any mode of any other panel are then derived. The 
general governing equations are obtained; these equations consist 
principally of three components: 
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-> panel modal structural impedances 


-> generalized forces due to the motion of a given 

structural mode of a given panel on another structural 
mode of another panel (including but not limited to the 
panel itself) 

-> a forcing term integrated over the panel mode shape 
which drives the system. 

The complexity of the general solution is discussed, along with the 
difficulties of obtaining a complete numerical and exact solution. A 
structural branch analysis is derived which includes the effects of 
higher structural modes. Clamped panels appear to be most stiff, 
whereas the flat rigid panel shape used in the models of Chapter 2 
and 3 appears the least stiff. A panel modal study shows that 
while it is theoretically possible to use higher structural modes to 
advantage in ART tuning real panels, these modes are often 
occurring at higher frequencies well out of range of practical use. 
Given the complexity resulting from using general panel mode 
shapes, the simple flat and rigid plate impedance model of Chapter 
2 and 3 yields conservative engineering results. 

Section 6.5: Exnerimenta l Confirmation of the AR1 
Concept 

A series of experiments confirmed the potential value of 
Alternate Resonance Tuning and displayed the relative ease with 
which the phenomenon can be obtained. A general experimental 
setup was constructed to investigate the ART tuning method for a 
number of practical geometries. Single panel arrays with as many 
as four panels and three ART frequencies show noise reductions of 
6 to 10 dB over the ANNR calculations. Measured phase 
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differences between passively driven audio speakers indicated 
that an approximately out of phase situation could be obtained for 
much of the frequency range between panel natural frequencies. 

Double wall experiments showed the possibility of using ART 
tuning in a multi-walled structure in numerous configurations, 
including nonuniform panel sizes. Real panel experiments 
indicated that the ART effect was not limited to simple devices 
such as passively driven speakers, but that legitimate panel 
sections could show ART noise reduction possibilities. Addition of 
a point mass to a Bakelite panel section “tricked” a plate into an 
apparent higher structural mode behavior which indicated ART 
cancellation potential of over 40 dB. 
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Section 6.6: Future Effort in Alternate Resonance Tuning 


Exploration of the ART concept continues at Duke University. 
A panel/frame method is currently being investigated for the 
solution of spatially periodic wall structures such as those found in 
aircraft^ 4 . A two-scale perturbation expansion method is used 
which is based on the scale separation between panel size and 
acoustic wavelength, and an averaging procedure is applied over a 
rapidly varying panel scale. A composite solution is formed which 
contains both a smoothed global solution and a periodic local 
solution that describes the details of the panel motions. Panel 
subsystems appear in an averaged sense in the global solution 
through transfer functions, obtained from the local solution for 
distributed frame loadings and acoustic boundary conditions. This 
panel/frame methodology holds great promise as a means to 
analyze complicated structures such as aircraft panel/frame 
assemblies, which appear “smeared” in a physical sense. 

Other more straight forward engineering analyses such as 
those performed in this research effort might be applied to 
problems to further understanding of the ART concept and its 
application to engineering structures. Figure 6.6.1 shows a 
schematic representation of a simple model for a multi-paneled 
wall structure such as those found in typical aircraft applications. 
This model considers a plate-like outer skin attached to a vibrating 
frame member, an inner rigid trim liner (also attached to the 
frame member), and a rigid plate hung on a suspension (much like 
an audio speaker) located between the two walls. The trim liner or 
suspended plate appears to be a practical way to implement ART 
tuning, and since ART appears to be most effect at low frequencies, 
the problem could focus on a two panel model to reduce 
transmitted noise at one frequency only. External panel tuning on 
an aircraft fuselage would be impractical; as mentioned in Chapter 
1, aircraft interior panels have been altered successfully in terms 
of mass and stiffness. Previous theoretical investigation of multi- 
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walled structures has encountered some difficulty in achieving the 
ART effect 35 , and more investigation may reveal new and 
innovative ways in which the ART concept might be used to help 
reduce noise transmission through paneled structures. 



Figure 6.6.1: Multi-walled model for future investigation. 
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